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ructed  waveform 


*  A  process  for  identifying  and  determining  the  structure,  depth 
and  relative  position  of  an  underground  tunnel  from  its  backseat tered 
response  is  presented.  Application  to  theoretical  and  measured 
tunnel  echoes  Is  included  with  very  satisfactory  results.  Basically, 
the  process  involves  the  Isolation  of  a  clutter  free  tunnel  response 
by  effectively  removing  all  undesired  poles  and  clutter  or  noise 
from  the  original  (unprocessed)  echo.  The  results  can  be  evaluated  — 
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Grey-Level  Mapping 
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by  Grey-level  Mapping  a  traverse  of  responses  over  the  tunnel.  The 
employed  s ignal-process ing  techniques  are  very  general  and  can  be 
>ised  for  detection  of  other  subsurface  targets. 
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CHAPTER  I 
INTRODUCTION 


A.  Background  and  Related  Research 

During  the  past  decade,  efforts  at  the  E lectroSclence  Laboratory 
have  been  concentrated  on  the  detection  and  Identification  of  various 
subsurface  targets  such  as  tunnels,  mine-like  targets  or  different 
kinds  of  pipes.  In  1965  Kennaugh  and  Moffatt  (1,21  performed  a 
pioneering  work  on  their  first  attempt  to  characterize  backscattered 
transient  responses  (echoes).  Young  and  Caldecott  13,4,51  later 
did  an  extensive  study  on  the  detection  of  pipes.  Recently,  Chan 
and  Peters  1 5-8 1  have  concentrated  on  the  detection  and  identification 
of  mine  targets.  Davis  and  Peters  (9)  give  a  short  summary  on  the 
tunnel  detection  progress  made  thus  far  at  the  E lectroSclence  Labora¬ 
tory.  The  above  contributed  greatly  to  the  advances  presented 
here. 

In  this  thesis  a  new  method  for  detecting  and  determining 
the  structure  (i.e.,  Ident ifying) ,  depth  and  relative  position  of 
underground  tunnels  is  presented.  It  was  applied  to  measured  tunnel 
backscattered  waveforms  obtained  at  a  tunnel  rite  located  at  Gold 
Hill,  Colorado.  The  results  as  compared  to  the  actual  tunnel  depth 
and  structure  were  reasonably  accurate. 

The  identification  process  is  based  on  the  assumption  that 
the  transient  response  of  the  tunnel  or  generally  any  target  can 
be  uniquely  character ized  by  a  set  of  complex  natural  resonances 
(poles)  (5,7,10.11).  We  excite  these  resonances  using  a  Video-Pulse 
Radar  (1?).  This  radar  consists  of  a  pulse  generator,  a  pair  of 
horizontal  crossed  dipole  antennas  (one  for  transmitting  and  the 
other  for  receiving,  thus  isolation  is  maintained  between  the  trans¬ 
mitted  and  received  signals),  and  a  samplinq  oscilloscope  is  used 
as  a  receiver.  The  generator  produces  a  narrow  pulse  of  broad  spec¬ 
trum  ranging  from  the  pulse  repetition  rate  of  the  pulser  to  the 
lower  microwave  region.  Although  the  transmitted  pulse  has  a  broad 
spectrin,  the  dipole  antenna  acts  as  a  principle  filter.  Therefore, 
the  natural  response  (or  simply  response)  of  the  target  must  contain 
its  resonances  close  to  the  dipole  resonance  for  possible  target 
detection  and  identification.  Tribuzl  and  Wald  (13,141  give  a  good 
description  on  the  development  of  the  various  dipole  antennas  that 
were  used  in  the  underground  radar  system.  It  was  found  that  the 
dipole  antenna  can  be  characterized  by  a  single  complex  conjugate 
pole  pair  (or  simply  pole  pair)  of  large  real  part,  independent 
of  antenna  position  (5-7). 
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The  oscilloscope  samples  the  received  echo  (of  some  finite 
time  window)  at  256  points  and  records  it  in  a  form  compatible  for 
computer  processing  1?  . 

B.  The  Tunnel  Response 

figure  1(a)  illustrates  a  simple  model  of  the  tunnel  response 
when  a  pulse  is  incident  on  its  surface.  The  incident  pulse  bounces 
back  and  forth  in  the  tunnel.  Every  time  the  bouncing  pulse  is 
at  the  upper  interface  of  the  tunnel  a  pulse  is  also  transmitted 
and  propagates  toward  the  antenna.  This  creates  a  multiple  lobe 
signal  structure  at  the  receiving  antenna  input.  Therefore,  (within 
our  system's  bandwidth)  the  tunnel's  (natural)  response  can  be 
character i/ed  in  the  cample*  frequency  plane  by  a  single  complex 
conjugate  pole  pair,  independent  of  antenna  position.  Figure  1(b) 
presents  a  transmission  line  model  of  the  unified  antenna-tunnel 
structure.  This  model  will  be  discussed  extensively  In  Chapter 
IV.  A  worst  case  study,  in  which  the  antenna  is  not  matched  to 
the  transmission  line  impedance  will  be  pursued  there.  In  practice, 
when  recording  actual  echoes  our  dipole  antenna  was  closely  matched 
to  the  ground  impedance,  thus  avoiding  multiple  reflections  from 
the  target  to  facilitate  the  processing  of  the  echoes,  especially 
when  dealing  with  shallow  targets.  The  results  of  this  theoretical 
study  will  prove  to  be  very  helpful  in  guiding  tunnel  identification 
process. 

To  avoid  confusion  in  later  discussions,  a  point  of  distinction 
is  due  here.  It  concerns  the  classification  of  the  different  wave¬ 
forms  to  be  encountered,  and  will  be  used  consistently  henceforth. 

When  we  refer  to  the  (backscattered)  original  waveform,  or  simply 
echo  we  mean  the  response  received  by  the  video-pulse  radar.  It  has 
not  been  processed  yet  by  any  identification  or  detection  scheme, 
except  some  analog  processing  accomplished  in  the  radar  itself  for 
supression  of  interference,  clutter  and  noise.  By  tunnel  or  target 
response  we  refer  to  the  natural  response  of  the  tunnel  or  target, 
respectively.  For  the  tunnel  this  is  just  an  exponentially  decaying 
s inusoid. 

C .  The  Identification  Process 

The  identification  process  is  concerned  with  the  calculation 
of  the  tunnel  resonance  and  the  arrival  time  of  the  tunnel  response. 
The  resonance  determines  the  tunnel's  height  and  the  arrival  time 
indicates  its  depth.  The  problem  associated  with  the  analysis 
is  that  the  received  backscattered  response  is  not  only  characterized 
by  the  tunnel  pole  pair  but  it  also  contains  the  antenna  pole  pair, 
possibly  other  false  target  poles,  and  significant  portions  of  clutter 
and  noise.  Clutter  occurrence  is  mainly  in  the  early  portion  of 
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the  received  echo  and  causes  difficulty  in  the  determination  of 
the  tunnel  response  arrival  time.  Major  efforts  are  attempted  in 
this  study  to  overcome  this  problem. 

In  processing  the  original  waveform,  there  are  several  distinct 
steps  and  the  curves  presented  are  labeled  accordingly.  These 
include: 


1.  Pole  fxtraction  Process 

This  is  in  reality  a  digital  or  time  domain  filtering  process 
where  a  natural  resonance  or  a  complex  conjugate  pole  pair  and  its 
associated  residue  are  removed  from  the  waveform.  The  reader  is 
cautioned  that  this  is  not  done  by  evaluating  the  residue  and  then 
subtracting  it  from  the  waveform,  but  instead  it  is  done  in  the 
time  domain  via  a  difference  equation  approach.  This  process  could 
indeed  be  referred  to  as  time  domain  filtering  and  indeed  for  the 
sake  of  clarity  and  simplicity  we  will  refer  to  it  as  such  on  the 
waveforms  presented.  The  approach  has  an  advantage  in  that  the 
natural  resonances  that  are  energised  at  different  times  can  be 
accounted  for. 

2.  Correction  Process 

In  the  process  of  extracting  a  natural  resonance,  the  remainder 
of  the  waveform  is  modified  by  this  filtering  process.  This  Step 
simply  takes  out  the  distortions  introduced  by  the  filtering  process. 

3.  Reconstruction  Process 

The  late  time  portion  of  the  signals  can  be  used  to  generate 
the  early  time  signal  of  a  given  pole  pair.  This  makes  it  possible 
to  reduce  the  effects  of  clutter  in  the  early  time  region  of  the  signal 
and  to  better  evaluate  the  initial  part  of  the  signal  reflected 
from  the  target. 

The  tunnel  identification  starts  by  first  recording  a  set 
of  backscattered  waveforms  from  a  pass  over  the  tunnel.  The  recorded 
echoes  are  first  processed  for  determining  the  poles  of  each  particu¬ 
lar  echo.  These  poles  will  include  the  antenna  and  tunnel  pole 
pairs,  and  possibly  other  false  target  resonances.  Prony's  classical 
method  (7,15-?01  was  used  for  calculating  the  poles.  This  was  done 
as  outlined  in  f7j,  by  applying  the  method  using  various  waveform 
time  windows  and  intervals  between  the  samples  to  be  used  for  con¬ 
structing  Prony’s  difference  equations.  The  selected  time  windows 
are  characterized  by  their  starting  time  and  length  (minimum  length 
is  dictated  by  the  manber  of  poles  desired).  According  to  the  cal¬ 
culated  poles  from  each  parameter  set  (window  starting  point,  window 
length,  manber  of  poles,  and  interval  between  samples)  a  theoretical 
waveform  is  constructed,  and  a  square  error  is  determined  between 
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In*  theoretical  and  measured  waveforms.  Finally,  the  poles  calculated 
by  the  parameter  set  giving  the  smallest  error  are  selected.  The 
program  “Singularity  Expansion  by  Prony“  as  given  by  Chan  [7]  was 
used.  Othei  techniques  such  as  the  Eigenvalue  method  [71  or  Contour 
Integration  |? 1 1  could  have  been  used  for  f indinq  the  poles. 

After  the  pole  calculation  we  strive  to  isolate  the  resoonse 
from  the  tunnel  alone.  This  is  the  most  important  step  of  our 
identification  process.  The  accuracy  of  our  results  as  compared 
to  other  attempts  by  Stapp  l??|  and  GEO-CENTERS,  INC.  [23 )  is  derived 
from  this  process.  It  is  accomplished  through  the  Pole  Extraction 
Process  and  is  extensively  discussed.  During  this  process,  all 
the  poles  not  associated  with  the  tunnel  are  removed  as  calculated 
by  Prony’s  method.  Thus,  ideally,  we  are  left  with  the  tunnel 
response  alone. 

Following  the  Pole  Extraction  Process  {anj  correction),  the 
Reconstruction  Process  is  used  for  combating  the  clutter  problem. 

It  reconstructs  the  early  portion  of  the  tunnel  response  based  on 
a  predicting  window  of  its  late  response.  This  process  contributes 
to  the  estimation  of  the  tunnel  response  arrival  time. 

Finally,  the  identification  process  can  be  highlighted  by 
mapping  the  reconstructed  set  of  tunnel  responses,  using  a  mapping 
technique  extensively  discussed  by  Stapp  [ 22 1 .  Such  a  map  can  indi¬ 
cate  the  relative  tunnel  position. 

0.  Structure  of  the  Thesis 

The  structure  of  this  thesis  is  as  follows: 

In  Chapter  II  the  Pole  Extraction  Process  is  discussed.  It 
presents  the  derivation  of  the  difference  equation  used  for  the 
extraction  of  undesired  poles.  Also,  various  character ist ics 
of  this  process  are  extensively  analyzed  and  criticized. 

In  Chapter  III  we  present  the  Reconstruct  ion  Process.  The 
difference  equation  used  for  accomplishing  this  process  is  derived, 
analyzed,  and  criticized.  Examples  on  the  performance  of  the  process 
with  measured  waveforms  are  given. 

In  Chapter  JV  a  simple  transmission  line  model  of  the  Radar- 
Tunnel  structure  is  analyzed  and  tested.  Application  of  the  Pole 
Extraction  Process  and  Reconstruction  Process  Is  given  here  for 
the  identification  of  the  theoretical  target. 

In  Chapter  V  measured  echoes  from  a  tunnel  are  encountered. 

The  effectiveness  of  the  Pole  Extraction  and  Reconstruction 
Processes  is  strongly  indicated  here  as  compared  to  other  attempts. 


chapter  II 

THE  POLE  EXTRACTION  PROCESS 


A.  Objectives 

The  objectives  of  this  chapter  are  the  following: 

1.  To  give  the  derivation  and  performance  of  the  difference 
»*quation  used  for  extraction  of  all  undesired  poles  from  the  raw 
recorded  waveforms  of  the  tunnel  echo. 

?.  To  discuss  the  effects  of  the  process  on  the  waveform 
associated  with  the  remaining  waveform  poles  and  the  importance 
of  the  sampling  interval  used  in  the  process. 

3.  To  derive  a  method  for  correcting  the  various  distortions 
occuring  on  the  waveform  associated  with  the  remaining  poles  due 
to  the  Pole  Extraction  Process. 

B .  Derivation  of  a  Difference  Equation  for  Extracting 

fine  Complex  Conjugate  PoV  Pair 

The  Pole  Extraction  Process  is  in  effect  a  filtering  process 
to  be  used  for  real  time  calcul at  ions .  It  is  accomplished  by  apply¬ 
ing  a  difference  equation  to  the  recorded  waveforms.  As  compared 
to  classical  filtering  it  works  in  the  complex  plane  (Figure  2) 
for  extracting  (removing)  particular  poles  from  the  recorded  echoes 
as  they  are  calculated  by  Prony's  method.  Its  main  advantage  is 
simplicity  and  speed.  It  avoids  convolution  or  frequency  spectrum 
calculations  as  would  be  required  with  classical  methods.  Furthermore, 
it  concentrates  only  on  the  particular  poles  to  be  extracted. 

In  order  to  understand  the  principle  of  the  process  let  us 
asstene  an  original  waveform  whose  spectrum  is  F(w)  and  is  character¬ 
ized  by  a  set  of  poles.  We  wish  to  extract  a  set  of  poles  charac¬ 
terized  by  a  function  Fj(w).  The  operation  of  the  pole  extraction 
process  is  as  follows: 

F?<»>  ■  r^T  •  (') 

F _ ( u)  is  the  spectrin  of  the  resultant  waveform.  Since  we  divided 
by  F.(io),  then  F~(u»)  has  all  the  poles  of  F(u>)  except  the  ones  of 
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Figurp  2.  Complex  plane. 


We  proceed  now  to  generate  a  difference  equation  for  imple¬ 
menting  a  process  **qu»valent  to  Equation  (1).  The  difference  equation 
for  extracting  one  complex  pole  pair  (2  poles)  Is  first  derived 
and  then  generalized  for  extracting  several  poles.  Since  we  are 
interested  in  deriving  a  process  to  be  applied  in  the  time  domain  and 
thus  avoid  spectrum  ca lcul at  ions,  all  manipulations  will  be  worked 
out  in  the  ^-domain.  £ -transform  representations  are  then  easily 
transformed  to  discrete  time  domain  |?5,?6).  Furthermore,  we  will 
only  stress  cases  dealing  with  complex  conjugate  pole  pairs  s^ce 
the  antenna  and  all  underground  targets  are  characterized  by  Such 
poles. 


0 

0 
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We  are  referring  here  to  one-sided  £  -transform  defined  as 

*• 

R(z)  •  r  r(nT)z'n,  where  R(z)  is  the  ^-transform  of  r(nT). 


]i»  t>) lowing  equivalences  between  Laplace  transforms  (contln 
uous  time)  and  ? -transforms  (discrete  time)  should  always  be  at 
hand: 


Sj  •  ot  ♦  J?»f  j 

\  *  J*, 

,  V 
*1  ’  e 

V 

/p  *  *  e  cos(Mf,T) 

'|T 

/y  *  !■(/,)  •  e  sin(?if.T) 

»  1 

°iT 

l*|l  "  e 


where  the  following  symbolic  notation  has  been  used: 

s  *  Laplace  transform  operator 

i  •  £  -transform  operator 

s4  «  pole  in  the  s-dcmain 

*  pole  in  the  /-domain,  equ i valent  to  s^ 

cij  ■  real  part  of  pole  Sj 

f  j  ■  frequency  of  pole  Sj 

/p  •  real  part  of  pole  zf 

/M  »  imaqinary  part  of  pole  z. 

1  1 
T  *  samplfnq  interval  of  the  process 

|z j |  *  maqnitude  of  / . 

/J  *  conjugate  of  /i 

s|  »  conjugate  of  s^. 


(2) 

(3) 

(4) 

(5) 

(6) 
(7) 
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The  above  not  jt  ion  will  be  consistently  used  in  all  subsequent 
references. 


let  ii*.  now  represent  our  original  waveform,  r(t),  in  the 
£  -  (lima  i  n  as  follows: 


R(*> 


_"UL_ 

yr.  w,  -r*: 


(8) 


where 
is  the 


(l-r'^Mi-r'zpou) 

R(z)  is  the  ^-transform  of  r(t)  (r(t)  >  (Zj.zj] 

■  cample*  pole  pair  to  be  extracted.  The  rest  of  the  denomin¬ 


ator  of  R(z),  D(z),  contains  the  other  poles  of  R ( z )  (tunnel  pole 
pair)  to  remain  after  the  extraction  of  (Zj.Zj). 


Multiplying  both^sides  of  Equation  (8)  by  the  representation 
of  the  pole  pair  (z,,Zj),  we  obtain 


y*>  *  Sffi  ■ 


(9) 


Rp(/)  is  the  ^-transform  of  the  desired  resultant  waveform.  Since 
we  have  mult ipl ied  R ( z )  by  the  zero. 


Rj (/)  * 

{ l-z'^^d-z 

the 

pole  pair 

of  R ( z )  at  (z 

contains  only 

the  poles  of 

(9) 

gives 

RpU)  • 

(l-z  Z  j -z  'z 

or 

Vz)  • 

( l-?Re(z1)z‘1 

(10) 


■ 2 .  •, 


,2.-2, 


(ID 


The  above  equation  can  now  be  easily  transformed  to  discrete  time 
domain,  to  obtain 


rp(nTe)-r(nTe)-?Re(z1)r(nTe-Te)H/1l?r(nTe-2Te) 


(12) 


Its  equivalent  form  In  continuous  time  Is  derived  In  Appendix  A. 
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where  I,  is  the, samp  I ing  interval  during  the  pole  extraction  process 
and  rp(t)  -  l«pU)]. 


The  sampling  interval  T^  is 


a  multiple  of  Tg,  i .e. , 


Ne  TU 


where  TR  is 
so  that  it 


the  basic  sampling  Interval  of  r(t). 
satisfies  Shannon's  sampling  theorem. 


N  must  be  chosen 
This  impl ies  that 


re  i  7T 


or  N 


(13) 


Equation  (1?)  as  it  stands  can  only  generate  points  which 
are  a  multiple  of  T  .  In  order  to  generate  all  the  points  within 
the  T  interval,  spaced  at  the  basic  interval  TR,  Equation  (12) 
can  be  modified  as  follows: 


VnVkV’r(nykV*?Re{*I)r{nVkVTe)  + 


I*/  r(nTe*kTB-27e) 


k*0.1 . Ne-1 

n*0, 1 ,2, . . . 


(14) 


The  above  equation  indicates  the  time  domain  operation  required  for 
extracting  the  complex  conjugate  pole  pair  (Zj,z.).  r  (t)  is  the 
residual  or  "filtered"  waveform  after  (*.,?.)  hai  beenPextracted. 
The  choice  of  the  proper  value  of  T  in  Equation  (14)  is  extremely 
important  and  the  implications  of  this  choice  will  be  discussed 
later  in  this  chapter. 

It  is  essential  to  observe  that  the  extraction  of  a  pole  pair 
requires  only  the  knowledge  of  the  pole  pair  itself  and  not  any 
information  about  its  residue.  This  is  of  great  advantage  since 
the  values  of  the  residues  are,  of  course,  excitation  dependent. 
From  Equation  (14)  it  is  noted  that  the  calculation  of  one  point 
of  r  (t)  requires  three  points  of  the  discrete  original  waveform. 
The  present  point  and  two  previous  points  spaced  at  intervals  of 
T  (or  NTg).  This  indicates  that  the  first  2N  points  of  the  fil¬ 
tered  waveform,  r  (t),  cannat  be  evaluated.  Ih  our  computer  imple¬ 
mentation  of  Equation  (14)  these  points  were  conveniently  set  equal 
to  zero. 
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The  method  for  implement ing  Equation  (14)  by  a  digital  computer 
is  shown  in  Figure  3.  As  seen  Equation  (14)  can  be  very  easily 
programm'd  by  a  field  microprocessor  for  real  time  calculations 
hy  simple  appropriate  shifts  and  additions. 


|  I  SAMPLES  OF  THE  ORIGINAL  WAVEFORM 

t 

•  :  SAMPLES  OF  THE  FILTEREO  WAVEFORM 
I 

Figure  3.  Computer  implementation  of  Equation  (14). 

An  application  of  Equation  (14)  on  a  theoretical  waveform 
is  shown  in  Figure  4.  The  original  waveform  (dotted  line)  is  com¬ 
pose^  of  two  exponentially  decaying  sinusoid^  with  polgs  at 
s.,Sj*-3.  Meqanepers/sec  ♦  j?«x15  Megarad/sec  and  s->,S2*-6.  nega- 
n^pecs/sec  ♦  J2xx20  Megarad/sec  and  residues  at  2.0  ♦  jO.O  for 
both  pole  pairs.  Its  representat ion  in  time  domain  is 

-6xl06nT  -3x106nT_  - 

r(nTg)  »  e  Bcos(2’ixl5xl0onTB)*e  Bcos(2"x20xl0t>nTB) . 

It  consists  of  256  points  with  T--200/255  nsec.  The  solid  line 
is  the  result  when  the  pole  at  l“  NH*  is  extracted.  As  seen,  it 
is  a  decaying  sinusoid  corresponding  to  the  s?  pole  pair  (20  HHz). 

The  pole  extraction  interval  is  Indicated  to  6e  10T„(N  ■ 10 ) .  Observe 
the  early  portion  of  the  filtered  waveform  which  is  set  equal  to 
zero  as  was  discussed  previously.  This  time  Interval  is  measured 
to  be  20Tg,  corresponding  to  15.7  nsec. 


♦  Indicates  the  presence  of  complex  conjugate  poles. 
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Figure  4.  Example  of  the  pole  extraction  process 
according  to  Equation  (14). 


I  ’  mj.mI  here  that  all  subsequent  waveforms,  measured  or 
theoretical,  will  consist  of  ?56  time  points. 

A  situation  of  importance  occurs  when  the  signal  associated 
with  one  of  the  pole  pairs  is  time  delayed.  This  is  usually  the 
case  for  echoes  from  deep  tunnels.  Then,  the  arrival  time  of  the 
tunnel  response  occurs  after  various  forms  of  clutter,  including 
direct  coupling  between  transmit  and  receive  antennas.  In  this 
case  when  the  antenna  pole  pair  is  extracted,  there  will  be  an  error 
region  of  time1  length  ?N  Tg,  corresponding  to  the  initial  portion 
of  the  delayed  part  of  the  echo  from  the  tunnel.  An  example  of 
such  a  situation  is  shown  in  Figure  5. 

C .  General  i/at  ion  of  the  Pole  Extraction  Process 
to  Several  Poles 


The  application  of  Equation  (14)  can  extract  only  one  complex 
conjugate  pole  pair  at  a  time.  A  similar  arrrrach  to  the  derivation 
of  Equation  (14)  can  be  used  for  deriving  difference  equations  to 
extract  concurrently  several  poles.  The  resulting  equations  will 
just  be  stated. 

The  difference  equation  for  concurrent  extraction  of  two  complex 
conjugate  pole  pairs  (4  poles)  or  less  can  be  shown  to  be  given 
by 


r 

P 


(nVkTB).r(nTeHTB).c,r(nTe+kTB-Te)4c2r(nTe*kTe-?Te) 


-c3r(nTe*kTB-3Te)*c4r(nTe»kTB-4Te) ;  k«0, 1 , . . . ,Ne-l 

n»0,l,2,... 

(15) 

where  r  (t),  r(t),  Ne,  T  and  T_  are  as  defined  previously,  and 
the  coefficients  are  calculate<rto  be: 

Cj  «  2|Re(/,)*ReU?)l 
c?  *  U)!?H*2l?*4Re(*1)Re(/?) 
c,  »  2l!/1|?Re(23M/?|?Re(*1)l 
C,  •  |i,l*l»2l*  . 

I 
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*  *  ~ 

and  (/-,./-)  -*re  the  j -domain  representat ion  of  the  complex 
conjugate  pole' pairs  to  be  extracted.  Since  four  poles  are  extracted, 
it  is  seen  from  Equation  (15)  that  the  calculation  of  a  single  point 
of  the  filtered  waveform  requires  four  previous  points  of  the  original 
waveform  and  its  present  point.  Therefore,  the  filtered  waveform 
cannot  be  calculated  in  an  initial  time  wondow  eaual  tc  4T  by  this 
difference  equation  when  two  complex  conjugate  pole  pairs  are  ex¬ 
tracted  concurrently. 

The  difference  equation  for  extracting  three  complex  conjugate 
pole  pairs  (6  poles)  or  less  is 


VnVk  V*r(nVkV"cir{nVkTB'V*c?r(nVkTB'?V 
-c3r<nVkV3V*c4r<nVkV4V 
'c5rfnVkV5V*c6r(nVkV6V;  k’0,1 . V1 

n*0,l,2,... 

(16) 

The  coefficients,  c^,  of  the  above  equation  are  given  as  follows: 
c,  *  2[Re(;1)*Re(/2)tRe)z3j] 

c?  *  !z| S?^ 1^2 1?^  i?^4Re(z7 )Re(z?)^4Re(i2)Re(z3)^4Re(z1)Re(23) 

c3  ’  2[]*i  l?R*(V*3>*l*2l?Re(V*3)*|*3l?Re(  VV* 
4Re(z1)ReU?)ReU;J)j 

c4  ■  l«il?l^|ZM*1|2^3l2*U2l?l?3l?^|2,|2Re(z2)Re(23)* 
4U?|?Re(z,)Re(z3)M|z3|?Re(z1)Re(z2) 

'S  •  I7 |Z3 12 1*, |*2 12 1*3 |?R*(i,)] 

c6  ■  U|l?U/l*3l? 
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where  (zj.zj),  (zp.z^)  and  ( Z3 . Z3 )  are  the  pole  pairs  to  be  extracted. 
It  is  seen,  that  the  coefficients  of  the  difference  equation  increase 
to  the  ntanber  of  poles  extracted  concurrently.  In  addition,  the 
nunber  of  previous  points  required  for  calculation  of  a  single  point 
is  also  equal  to  the  nixnber  of  extracted  poles.  This  necessitates 
the  expansion  of  the  error  region  of  the  filtered  waveform,  equal 
to  ml  when  m  complex  poles  are  extracted  concurrently,  larger  time 
window^  of  the  hriginal  waveform  are  then  required  or  smaller  T 
if  we  are  to  use  r  ( t )  in  subsequent  processings.  This  problem e 
'  an  be  alleviated  by  use  of  toe  Reconstruct  ion  Process  discussed 
in  the  next  chapter. 

The  computer  subroutine  given  in  Appendix  C  implements  Equation 
(16).  An  application  of  Equation  (16)  is  presented  in  Figure  6. 

The  results  were  produced  by  the  main  program  given  in  Appendix 
B  which  also  uses  the  pole  extraction  subroutine.  The  dotted  line 
in  Figure  6  is  a  theoretical  wavefonp  composed  of  three  complex 
conjugate  pole  pjjirs  located  at  S|,sj*-6.  ^ganepers/sec  ♦  j?nx50. 
»**qa[ad/sec,  s^s?*-6.  teganepers/sec  ♦  j?vx40-  Megarad/sec* 
s3*s1*'6.  Meganepers/sec  ♦  j?»x?0.  Megarad/sec ,  with  respective 
residues  at  1 . ♦JO . .  1 . ♦ JO  and#?.*jO.  The^solid  line  Is  the  result 
after  the  extraction  of  ( s 1 , s ^ )  and  (sp.s^).  It  is  jimply  an  ex¬ 
ponentially  decaying  sinusoid  corresponding  to  (53,53).  The  coef¬ 
ficients  are  calculated  to  be: 


c,  ■  2.  46/69! 
c?  *  3.476065 
C3  *  7.376516 
c4  *  .92741) 


for  Np*4.  Note  that  the  error  region  is  equal  to  4NeTg« 1 ? . 55  nsec. 

A  generalized  difference  equat ion  can  be  found  for  extracting 
any  njwber  of  poles,  m  .  The  $ -transform  representation  of  the 
filtered  waveform  after  the  extraction  of  m  poles  can  be  given 
as  e 


■  *u)  TT  d-*,*"1)  (i7) 

M 

where  z.  represents  the  extracted  pole.  If  z,  is  complex,  then 
the  product  must  also  include  its  complex  conjugate  when  working 


Figure  6.  Application  of  Equation  (16);  concurrent  extraction  of 
two  complex  conjugate  pole  pairs. 


with  real  waveforms.  Expanding  the  product  of  Equation  (17)  (24] 
and  transforming  to  the  time  domain  we  obtain  the  following  general 
pole  extraction  difference  equation: 

m 

V"V*V'r<"VkTB>‘<-'>?  I ,  ;t1r("VkW 

1^1  1 


e  6 


<•"  1,1.  kl/<"VkVi!T.> 


1,-1  t?«l  *1  -2 

Vl2 


,3  r° 


m 

e 

m 

e 

l 

l 

V1 

V1 

■«3 

1 

e 

(-D 


7  l  i  l  r(nT  ♦kT_-3T  ) 

i-  l.  l,  li  1  e  B  e# 


ri  •••*•  r(nVkV1Te) 

eg  •  M  e  o  e 

i ,  1  1  c 


rn 


ll,t?' 


>v, 


6 

k*0, 1 ,2, . . . ,N  - 1 

n ■0,1,2,...  (18)  + 


The  ? -transform  representat ion  of  this  equation  also  appears  in 
ChapTer  IV  of  Chan  ’7'. 


18 


I  iio  reader  can  easily  varity  that  Equation  (  I H )  can  be  readily  reduced 
to  equations  (14),  (IS)  or  (16). 


D.  effects  of  the  Pole  Extraction  Process 

or.  f^e  yiveforwi  Xssoc  i  ated  With  the 

Rem .1  TriT n^  Poles 

We  now  turn  our  attention  to  another  important  parameter 
ot  the  pole  extraction  process,  its  sampling  interval,  T  .  When 
Equation  (18)  is  used  to  extract  the  desired  poles  then  the  residues 
associated  with  the  remainder  of  the  waveform  are  changed  in  both 
magnitude  and  phase.  These  modifications  will  be  shown  to  be  a 
main  function  of  T  and  also  the  natural  resonances.  To  acquire 
a  feeling  for  thesA  effects  we  refgr  to  Figures  7  and  8.  The  original 
waveform  is  composed#of  poles  S] ,Sj*-6.  Meg a nepers/ sec  ♦  j^nxl5. 
M"garad/sec  and  sp,s^»-6.  Mpganepers/sec  ♦4j?nx?0.  Megarad/sec  with 
residues  of  in  both  figures  (si.Sj)  is  extracted.  The 

extraction  interval  of  Figure  7  is  5Tg(Tg*500/255  nsec)  and  that 
of  Figure  8  is  9Tg.#  The  filtered  waveform  is  characterized  by 
the  pole  pair  (St.s^),  but  its  magnitude  and  phase  is  different 
for  each  extraction' interval .  In  fact,  none  of  these  results  corres¬ 
pond  to  the  proper  residue  of  the  original  signal  associated  with 
the  pole  pair  (s^.s?)  (its  proper  form  is  also  shown  for  comparison). 
The  case  when  the  original  waveform  is  composed  of  poles  Sj,Sj*-6. 
^•ganepers/sec  ♦  j?  ix30.  Megarad/sec  and  s2,S9*-6.  Weganepers/sec 
♦j?-'x?0.  Megarad/sec  with  residues  of  ?.♦  j6.  is  shown  in  Figure  9. 

The  extraction  interval  is  5Tg  (same  as  in  Figure  7).  As  seen, 
the  residue  of  the  remaining  signal  does  not  have  the  same  amplitude 
ir  phase  is  that  of  F  iqur«>  \  The  measured  amplitude  and  phase 
lis^ortions  are  indicated  on  the  respective  figures. 

In  order  to  quant i t at i vely  present  the  effects  of  T  and  the 
extracted  pole  pairs  on  the  remaining  signal,  a  simple  2-Pole  pair 
(complex  conjugate)  waveform  of  the  form 


r(nTB)*?|A1 


OinTg  OphTg 

e  cos(u>jnTg'Kij )  ♦?  |  A^|  e  cosfu^nTg+a^) 


(19) 


is  first  studied.  We  will  then  generalize  our  results  to  include 
additional  poles.  The  ^  -transform  of  Equation  (19)  is 


RU) 


‘1 


A,  A? 

rr*  *  7T-T ;  * 


T7rrz1  I-*-1**  1-z‘^Zp  Uz‘Tz  ^ 


(20) 


Jo, 


where  Sj.s,  •  OjlJuj.  s2,s2*o2iJu?,  A^jA^e  and  A?*|A?|e 


39V110A 


Figure  7.  Example  on  the  amplitude  and  phase  effects  due  to  th 
extraction  process. 


pole  extraction  process. 


The  extraction  of  pole  pair  (*.,*.)  Involves  the  multiplication 
of  R(z)  by  Equation  (10).  Then,  the  transform  of  the  filtered  wavefo 
is  given  by 


Aj(l-»e(*1)r,*|*,|V2) 

^-rrv 

A2(l-?ReU,)r,*|i1|2*"2) 

l-Z^Z, 


Expanding  the  last  term  of  the  above  equation,  we  have 


Performing  long  divisions,  we  find  that 


*2r  i*2ro-*"*2) 


l*2rn-*  ’*?) 


The  expansion  of  the  t$1rd  term  of  Equation  (fl)  Is  similar  to  the 
above.  If  A,,  z,  and  z,  are  replaced  by  A,,  z,  and  z9,  respectively. 


After  substitution  of  the  above  expansions  in  Equation  (21), 
we  can  express  Rp(z)  as  follows: 


V2) 


A,Z.,*A,Z,\  ,  f A,( z,)^*A,( Z-)^ 

»t<A,).».(z,)  -  |z,|2  -2-*— f-*- 

„  'z2'  /  \  <*2^ 


-  1  •  •  I  ^  1  I  *  * 

-  *  A,z 1  ♦  A,?,  ♦  -  (A?z?*A?z?) 


l-*"z2  l-z"z2 


f—  *  «  ?  “ 

♦  -?Re(z j)  —7  ♦  Izjl  —4—  T“rr" 

I *2 '  1*^1  1 -z 

.  ,  V*  .  |f  l*  ■it'afl— » 

1  i«jp  1  iz2r 


The  first  and  second  term  of  Rp(z)  in  Equation  (25)  corresponds 
to  impulses  at  t«0  and  t»T  *N jL,  respectively.  As  was  Indicated 
previously,  these  points  cannotMbe  calculated  by  Equation  (12). 
Therefore,  they  will  be  discarded  in  the  rest  of  our  discussion. 

The  third  term  corresponds  to  the  original  pole  pair,  (Zj.Zp). 

This  is  the  desired  result.  The  fourth  and  fifth  terms  are ^distor¬ 
tion  terms  and  correspond  to  the  variations  on  the  magnitude  and 
phase  of  the  residue  for  the  signal  associated  with  the  remaining 
pole  pair.  Therefore,  R„( z )  can  effectivejy  be  written  as  the  sum 
of  the  original  signal  with  poles  at  (z-j.Zp)  and  the  distortion 
terms.  Note  that  these  terms  involve  z,  and  z,  which  are  functions 
of  T  .  >2 


In  order  to  understand  the  effects  of  the  distortion  terms 
let  us  proceed  by  first  inverting  them  into  time  domain.  Thus, 


Distortion  Tertns 


ReU,)  ^le  J2(nTe*kTB)  fJt2(nVkTB,4a2'Ve] 


-2|A?|  - e  "  'e 

i1?' 


.  -J[“?(",e*kT8)*',2-u2T,] 


.|S,|  !fi£  /°2V2l"I»*k7B)  (/“><"  VkV*V2V,] 

*J  :<^(nTeHV*V?u,2Te](  «k*®*  ^  •  •  •  • 

*  e  U  J  n-0,1,...  (26) 


Or 


Distortion  Tertns  * 


(o.-o?)T  °?(nT  ♦kT.)  r- 

-4|A?|e  ecos(«jTe)e  *  <^2(nVkV*VVeJ 


2(o.-0-)T  o.(nT  *kT  )  r  , 

♦2]A?|e  1  ee  1  e  cos!*y  "VkV*  V2U,2TeJ; 


k»0,1,...,Ne-1 
n»0, 1 , . . .  .  (27) 

For,our  waveforms,  typical  values  of  o.  and  o?  are  in  order  of 
-10°  Nepers/sec,  while  those  of  T  in  10*y  nsec.  Me  can  then 
use  the  approximation  e 

o  t  o  t 

e,e»e?esl.  (28) 


According  to  this  approximation,  r  (t)  (Equation  (25))  can  he  ex¬ 
pressed  as  follows:  9 


(°1'°?)TP  °?(nTPfkTR) 

r p(nVkT8)  8  2lA?!e  1  *  1 


cos[u2(nTe*kTB)^a2j 


(ciro?)T  o?(nT  ♦kTg)  r  *1 

-4|A?|e  e  i  cos(u>1Te)cos^i?(nTe*kTB)'«»?-<J2TeJ 


(o.-o«)T  u«(nT  ♦kTp)  <-  —i 

►?|A?|e  £  e  '  e  cos,u,2(nTe*kTB)^a2-?u.2TeJ  ; 

k*0, 1 , . . . ,Ne*l 
n»0,l,...  .  (?9) 


After  further  manipulation,  by  using  some  of  the  cosine  identities, 
we  arrive  at  the  important  result: 

o9("T  ♦kT„-TJ  - 

rp(nVlcV84lA?le  ‘  cosj^2(nTe+kTB)*a2-u>2TeJ  x 

o  ■  T  * —  -j 

e  cos(a2Te)-cos{u>,Te)J  ;  k=0,1 . Np-1 

n*0, 1 . 

(30) 

We  conclude  that  the  signal  associated  with  the  remaining  complex 
conjugate  pole  pair  will  be  shifted  in  phase  by 


-es.W2Te-«2NeTB-2»f2NeTB(rad) 


(31) 


or  equivalently  will  be  time  increased  by 
ft 

t  *  * — r-  •  HI-  *  T  , 

s  ?*1 ^  e  8  e 

It  will  also  be  amplitude  distorted  by  a  factor 
AF *?e  ,Te[cos(?*f?NeTB)-cos(?’'f1NeTB)]  . 


(3?) 


(33) 


?6 


We  observe  that  the  time  shift  does  not  depend  on  the  inter¬ 
acting  pole  pairs.  It  is  equal  to  the  extraction  interval.  Con¬ 
versely  the  amplitude  factor  does  depend  upon  the  frequencies  of 
the  pole  pairs  and  also  T  .  If  it  is  negative,  this  is  equivalent 
to  an  additional  phase  shift  of  A  «-*  or  a  time  Increase  of 
1  5 

t$  ■  yjr-  .  When  sampling  at  the  Nyquist  rate  (Equation  (13))  of 

f^  then  0  is  multiples  of  ».  The  amplitude  factor  (Equation  (33)) 
is  plotted  in  Figure  10  as  a  function  of  f-.  It  is  interesting 


T  «  200/255  NSEC 


NjtU'JtN 


Amplitude  response  of  the  pole  extraction 
process  (Equation  (33))  vs.  frequency 
(f^)  of  the  remaining  pole. 


to  note  that  the  amplitude  of  the  filtered  waveform  can  increase 
fourfold  when  ijpNeTg»?nx  and  ‘**]NeTg»(?n*l)w,  and  conversely.  But 


i.  .-n  in  T  I  (jure  10,  it  can  also  vanish  to  zero  when  w?N^Tg*2nn 
and  ?k».  Therefore,  we  must  always  select  a  sampling  interval 

which  does  not  bring  us  close  to  the  zero  crossings  of  the  ampli¬ 
tude  response.  In  addition,  when  N  is  small,  then  the  two  cosine 
terms  approach  1.  Thus,  the  .amplitude  factor  tends  to  zero.  This 
is  the  reason  for  performing  the  pole  extraction  process  at  a  samp¬ 
ling  interval  close  to  the  Nyquist  interval  of  the  extracted  pole 
pair  with  the  highest  frequency  (7). 

It  should  be  noted  that  Equation  (30)  is  not  valid  when  z,*z,. 
For  this  case  Equation  (20)  must  be  modified  to  treat  double  poles'; 

It  can  be  easily  shown  that  when  double  poles  are  encountered,  there 
is  no  distortion.  After  extraction  of  one  of  the  complex  conjugate 
pole  pairs,  the  signal  associated  with  the  remaining  single  pole 
pair  is  returned  in  its  original  form. 

Furthermore,  Equation  (33)  can  be  studied  for  the  case  when 
the  extracted  pole  pair  is  slightly  different  than  its  actual  value 
in  the  waveform.  This  is  generally  the  case  with  measured  waveforms. 
The  poles  calculated  by  Prony's  method  are  an  approximation  of  the 
actual  ones.  The  accuracy  of  the  results  depend  on  the  noise  and 
clutter  level  present  in  the  waveform  [7,18;.  From  Equation  (33) 
we  observe  that  if  the  deviation  in  frequency  is  small  there  is 
no  great  deter lorat ion  of  the  process.  The  cosine  terms  are  almost 
equal  and  the  amplitude  factor  becomes  negligible  (Equation  (30) 
goes  to  zero).  Any  deviation  in  the  real  part  of  the  pole  pair 
is  less  critical  (see  Equation  (28)).  If  it  is  large,  then  the 
desired  pole  pair  will  not  be  extracted  completely,  but  it  will 
be  reduced  in  magnitude,  as  in  the  case  of  frequency  deviation. 

When  the  measured  data  are  very  noisy  it  would  be  helpful  if  the 
poles  are  calculated  again  in  order  to  evaluate  the  pole  extraction 
process.  It  might  be  necessary  that  the  pole  extraction  process 
be  repeated  again  according  to  the  new  calculated  poles. 

In  order  to  evaluate  our  results  in  Equations  (32)  and  (33) 
we  can  compare  them  with  the  measured  ones  in  Fiqures  7,  8  and  9. 
Referring  to  Figure  7  according  to  Equations  (32)  and  (33)  we  haves 


ts»9.8  nsec,  AF*-.51 

or 

ts*38.8  nsec,  AF=.51  . 

Similarly,  for  Figures  8  and  9  we  calculate; 

^=42. 6  nsec,  AF*.92 
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*,.*9.8  nsec. 


AF*  1 . 14, 


respectively.  As  seen.  Equations  (32)  and  (33)  perfectly  predict 
the  t imo  shift  and  amplitude  distortions. 

E .  Correction  Process  Involving  Interaction  of 

Two  Complex  Conjugate  Pole  Pairs 

The  amplitude  and  phase  distortions  on  the  signal  associated 
with  the  remaining  pole  pair  due  to  the  pole  extraction  process 
is  an  undesirable  effect.  We  want  to  obtain  the  remaining  portion 
of  the  waveform  in  its  original  form.  Otherwise,  errors  will  occur 
in  the  tunnel  identification  process  to  follow.  It  will  be  seen 
in  Chapters  IV  and  V  that  if  the  resultant  target  response  is  not 
properly  time  placed,  then  errors  will  occur  when  determining  the 
target's  depth  and  position.  Furthermore,  amplitude  distortions 
could  obviously  diminish  the  original  tunnel  response. 

It  is  essential,  then,  to  correct  for  the  amplitude  and  phase 
distortions  after  the  pole  extraction  process  has  been  completed. 
This  is  easily  accomplished,  using  Equation  (30),  as  follows: 


rc<t)  *  rpU*ts)/AF.  (34) 

where  r ^ ( t )  is  the  corrected  waveform. 

The  correction  process  of  Equation  (34)  was  applied  to  the 
filtered  waveforms  in  Figures  7,  8  and  9.  The  result  (solid  line) 
after  correction  is  shown  in  Figures  11,  12  and  13,  respectively. 

For  comparison,  the  actual  form  of  the  signal  associated  with  the 
original  pole  pair  (s2,s2)  is  also  plotted.  As  seen  there  is  a  perfect 
match  of  t(je  corrected  and  original  signal  associated  with  pole 
pair  (s2,s?). 

The  correction  process  will  be  Incorporated  in  all  subsequent 
pole  extractions.  When  we  refer  to  filtered  waveforms  it  is  assumed 
that  correction  has  been  already  performed. 

F.  Correction  Process  Involving  Multiple 
(Complex  Conjugate)  ftoTe  Pairs 

In  the  previous  section  we  discussed  the  correction  process 
(Equation  (34))  when  two  pole  pairs  (complex  conjugate)  were  Involved. 
Obviously,  measured  waveforms  will  not  necessarily  consist  of  only 
two  pole  pairs.  Two  questions  can  then  be  raised:  How  is  the  re- 
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Figure  11.  Application  of  the  correction  process  to 
the  result  of  Figure  7. 


Figure  12.  Application  of  the  correction  process  to 
the  result  of  Figure  8. 


Figure  13.  Application  of  the  correction  process  to 
the  result  of  Figure  9. 


mainmg  waveform  effected  when  two  or  more  pole  pairs  are  extracted 
concurrently,  and,  how  do  we  account  for  the  case  of  a  remaining 
signal  associated  with  more  than  one  pole  pair? 

The  first  question  is  easily  answered  by  observing  that  each 
extracted  pole  pair  will  affect  the  signal  associated  with  each 
remaining  pole  independently  according  to  Equations  (3?)  and  (33). 
Therefore,  the  correction  process  is  performed  by  the  application 
of  Equation  (34)  for  each  of  the  extracted  pole  pairs  on  the  re¬ 
maining.  This  type  of  multiple  correction  is  included  in  the  pole 
extraction  subroutine  in  Appendix  C.  An  example  of  multiple  pole 
pair  extraction  and  correction  is  shown  in  Figure  14,  corresponding 
to  the  uncorrected  case  gf  Figure  6.  The  original  signal  associated 
with  the  pole  pair  ^.Sp)  is  also  plotted  for  comparison. 

The  second  question  can  be  difficult  since  the  correction 
process  requires  the  separation  of  the  signal  associated  with  the 
remaining  waveform  pole  pairs.  This  case  is  not  encountered  in 
tunnel  identification,  since  the  tunnel  can  be  characterued  by 
a  single  pole  pair.  A  method  for  performing  such  an  operation  is 
diagrammed  in  Figure  15.  We  have  an  original  waveform  consisting 
of  pole  pairs  A,  B,  C  and  0.  We  wish  to  extract  only  A  and  B. 

Pole  pairs  C  and  D  are  to  remain.  The  pole  extraction  process  takes 
two  directions  in  this  case.  In  one  direction  the  signal  of  pole 
pair  C  is  isolated  and  corrected  appropriately,  and  in  the  other 
the  signal  of  pole  pair  D  is  treated  likewise.  The  last  step  involves 
the  addition  of  the  corrected  signals  with  pole  pairs  C  and  D. 

This  approach  is  effective  on  a  theoretical  basis.  The  author  has 
not  exploited  the  method  when  dealing  with  measured  waveforms. 
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Figure  14.  Pole  extraction  and  correction  involving 
multiple  extracted  pole  pairs. 


Figure  IS.  Pole  extraction  and  correction  with  multiple 
extracted  and  remaining  complex  conjugate 
pole  pairs. 
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CHAPTER  III 

THE  RECONSTRUCTION  PROCESS 


A.  Object ives 

The  objectives  of  this  chapter  are  the  following: 

I.  To  discuss  the  problems  associated  with  clutter  and  clipping 
in  the  backseat tered  waveform  from  the  tunnel. 

?.  To  derive  and  discuss  the  proper  application,  and  advan¬ 
tages  of  the  Reconstruction  Process  used  for  overcoming  the  clutter 
problem.  This  involves  a  simple  difference1  equation  which  uses 
a  portion  of  the  late  tunnel  response  to  reconstruct  its  earlier 
time  portion.  A  criterion  for  selecting  this  late  portion  of  the 
tunnel  response  is  also  presented.  For  the  sake  of  completeness 
the  reconstruction  process  is  generalized  to  reconstruct  target 
responses  character i zed  with  more  than  one  pole  pair. 

3.  To  apply  the  reconstruct  ion  process  on  measured  tunnel 
responses.  This  is  accomplished  by  first  extracting  the  antenna 
pole  pair  and  possibly  other  false  target  resonances.  The  char¬ 
acter  of  the  residual  signal  after  the  pole  extraction  process  is 
classified  and  discussed  for  the  correct  and  effective  application 
of  the  reconstruction  process  (Figure  17).  It  will  be  indicated 
that  the  later  t tme  of  the  filtered  waveform  will  contain  the  tunnel 
response  exclusively,  except  for  clutter  and  noise.  Therefore, 
the  portion  of  this,  under  a  mean  square  error  criterion,  is  Chosen 
for  reconstructing  the  tunnel  response.  Thus,  we  obtain  a  "clutter 
and  noise  free"  tunnel  response  if  we  can  Indeed  find  such  a  window 
in  the  late  tunnel  response.  The  reconstructed  tunnel  response 
will  be  used  later  for  estimating  the  tunnel  depth. 

B.  The  Clutter  Problem  j 

In  the  introduction  it  was  noted  that  the  early  portion 
of  the  backscattered  echoes  were  dominated  by  clutter.  This  is  j 

mainly  due  to  ground  reflections.  Such  clutter  creates  high  voltage 
peaks  which  were  clipped  before  processing  by  the  recording  oscilo-  j 

scope.  Typical  tunnel  backseat tered  echoes  are  shown  in  Figure 
16.  These  were  obtained  from  the  tunnel  site  to  be  discussed  in 
Chapter  V.  Their  late  time  portion  was  found  by  Prony's  method 
to  contain  two  natural  resonances  or  pole  pairs  (complex  conjugate), 
one  corresponding  to  the  antenna  and  the  other  to  the  tunnel. 
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Appearance  of  clutter  in  the  early  portion  of  the  recorded 
waveform  disrupts  useful  information.  The  most  important  one  is 
the  arrival  time  of  the  tunnel  response.  Knowledge  of  this  is 
essential  since  it  will  enable  us  to  estimate  the  distance  of  the 
tunnel's  top  from  the  antenna.  Furthermore,  maps  constructed  from 
echoes  of  a  set  of  measurements  when  the  radar  is  plac^  above  the 
tunnel  become  difficult  to  read  and  the  region  corresponding  to 
the  tunnel  depth  cannot  be  identified  unless  some  previous  knowledge 
about  the  depth  is  given.  Basically,  clutter  is  an  undesired  effect 
that  contains  no  apparent  useful  data.  In  addition  it  disrupts 
the  good  data  to  follow. 

Stapp  12?)  and  a  group  at  GEO-CENTERS,  INC.  attempted  to  con¬ 
struct  maps  without  accounting  for  the  clutter  interference  in  the 
recorded  waveforms.  Their  maps  are  confused  and  do  not  clearly 
contain  the  distinct  tunnel  behavior  as  a  function  of  radar  position 
when  the  radar  data  are  used  to  develop  a  map.  It  is  quite  difficult 
to  obtain  precise  information  about  the  target  using  their 
maps. 


C.  Derivation  of  a  Difference  Equation  for 

Reconstruct Tnq  the  Early  Tunnel  Response 

A  method  is  introduced  here  to  overcome  the  clutter  problem. 

It  is  used  in  conjunction  with  the  pole  extraction  process  and  it 
will  be  referred  to  as  Reconstruction  Process.  It  is  applied  after 
ALL  undesired  poles  (mainly  the  antenna  pole  pair)  have  been  extracted 
from  the  original  waveform.  Thus,  a  "tunnel  response  alone"  is  ob¬ 
tained.  The  reconstruction  process  assimes  that  any  information 
contained  in  the  early  portion  of  the  recorded  waveform  cannot  be 
extracted.  Therefore,  it  is  discarded.  In  consequence  the  late 
tunnel  response  is  used  to  reconstruct  its  early  response  based 
on  the  knowledge  of  the  tunnel  pole  pair  as  determined  by  Prony's 
method.  This  is  accomplished  by  applying  a  difference  equation 
to  predict  early  points  from  later  ones.  The  derivation  of  the 
difference  equation  used  for  reconstruction  will  now  follow: 


Let  us  assume  an  original  echo  consisting  of  m  *2  complex 
poles.  The  two  correspond  to  the  tunnel  pole  pair  and  the  other 
m  are  undesired.  After  their  removal,  the  A  -transform  of  the 
filtered  waveform  (it  was  noted  at  the  end  of  the  last  chapter  that 
all  filtered  waveforms  are  assigned  to  be  corrected  according  to 
Equation  (3A).  The  reconstruction  process  has  no  usefulness  unless 
the  filtered  waveform  is  corrected  before  reconstruction)  can  be 
expressed  as 


w 


-1 


R  (i)  *  - 7T~ - y 3 

l-2Re(z,)*  '♦|z1!Zz  Z 


(35) 
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where  ( z  | .  2  ^ )  is  the  remaining  tunnel  pole  pair.  Multiplying  both 
sides  of  Equation  (35)  by  the  denominator  of  the  fraction  on  the 
right  hand  side  gives 

(l-?Re(21)2’14|z)|?z-?)Rc(i)  .  Va,*’1  •  (36) 

The  rjght  hand  side  of  Equation  (36)  is  a  first  degree  polynomial 
in  i~  .  Thus  it  can  be  expressed  as 

(l-?Re(71)/',4|z1|?z'2)Rc(i).A(;)  .  (37) 

Transforming  Equation  (37)  to  a  discrete  time  function,  we  obtain 


rc{nVkTB)  -  )rc(nTr4|lT8~Tr)4 
l*ll?rC<"VkV2Tr>  *  rA(nVkV 


k*0, 1 . Nr*l , 

n*0,l,...  (38) 


where  r,(nT  )  "  3  "'lA(*)l*  Note  that  Tr  again  is  not  the  basic 
sampling  interval  of  the  original  waveform,  T„.  It  is  a  multiple 

T  1  _  B 
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The  important  observation  to  be  made  in  Equation  (38)  is  that: 


rA(nTr)»0  for  n  _>  2 

Therefore, 


(40) 


rC(nTr*kTB)'?ReUl)rc(nTr4kTB‘Tr)4,Zll?rc(nVkTB‘2Tr)'0; 


"i?,  **0.1 . Np- 1  .  (41) 


Rearranging  the  above  equation,  gives 


rc<"V“V?Tr>- 


?Re(/,) 


1 


rc(nTr*kV; 


This  equation  can  predict  an  early  point  based  on  the  knowledge 
of  two  later  points  and  the  (tunnel)  pole  pair  associated  with  the 
portion  of  the  waveform  to  be  reconstructed.  The  equation  is  only 
good  when  applied  in  the  region  where  the  tunnel  pole  pair  is  cor¬ 
rectly  described.  Therefore,  it  is  true  for  the  region 


nTr  1  WeV2Tr  (43) 

where  t  represents  the  starting  time  of  the  portion  of  the  waveform 
associated  with  the  desired  tarqet,  i.e.,  the  tunnel  as  shall  be 
discussed  later  in  detail.  As  before  N  T„  is  the  extraction  interval 
used  when  extracting  the  m^  undesired  cSmplex  poles.  The  term 
m  NeTg  accounts  for  the  error  region  due  to  clutter  and  clipping. 

These  points  are  better  illustrated  in  Figure  17.  There,  the  antenna 
pole  pair  [-27.  Meganepers/sec  ♦  J2*x45.5  Megarad/sec]  of  the  echo  in 
Figure  16(b)  is  extracted.  The  respective  regions  are  indicated  on 
the  graph.  The  late  portion  of  the  filtered  waveform  contains  the 
scattered  field  from  the  target  (tunnel)  almost  exclusively.  There 
are,  of  course,  other  signals  that  could  be  present,  caused  by  clutter 
and  other  scatterers  if  their  resonances  are  within  the  radar  band¬ 
width.  In  this  case  the  poles  of  these  scatterers  (false  targets) 
must  be  removed  by  the  extraction  process  as  was  done  with  the  antenna 
pole  pair  for  obtaining  the  target  response  exclusively. 

Note  that  in  Equation  (42)  newly  generated  sample  points  can 
be  used  for  reconstructing  earlier  ones.  Therefore,  only  an  initial 
“errorless"  time  window  of  length  ?N  T.  is  needed  for  the  recon¬ 
struction  process  to  begin.  In  practice  (see  Figure  17)  we  actually 
never  encounter  an  "errorless"  time  window.  For  this  case  we  strive 
to  choose  the  window  that  best  describes  the  tunnel  pole  pair  calcu¬ 
lated  by  Prony's  method.  This  is  the  subject  of  the  next  section. 

The  implementation  of  Equation  (42)  is  similar  to  that  shown  in 
Figure  3,  except  that  reconstructed  samples  are  used  for  the  region 
out  of  the  chosen  initial  time  window. 

Furthermore,  It  is  essential  to  realize  that  Equation  (42) 
does  not  assume  the  presence  of  any  dc  term  in  the  filtered  waveform. 
It  is  imperative  then  that  any  dc  is  removed  before  the  process 
is  applied. 
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D. 


Error  Criterion  for  Selecting  the 
Base  T ime  Window 


The  problem  to  be  encountered  now  Is  the  selection  of  the 
initial  (base)  time  window  to  be  used  for  reconstruct  ion.  The 
criterion  for  such  a  decision  was  chosen  so  that  the  first  recon¬ 
structed  points  must  best  fit  the  filtered  waveform.  In  this  case 
the  reconstruct  ion  process  is  applied  by  using  different  base  windows. 
Actually,  the  portion  of  the  filtered  waveform  beyond  the  clutter 
region  is  scanned  by  shift inq  the  base  window  from  right  to  left 
by  Tg.  For  each  base  window  used,  an  error,  er,  is  calculated  as 


256-N 

V  ^ 
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er  * 


i«?S7-Ns-k 


rc(1TB)-rr(’V 


(44) 


where  r  (iTg)  is  the  filtered  waveform,  r  (iT«)  is  the  reconstructed 
waveform,  fr  is  the  reconstruction  starting  point,  counted  from 
the  end  of  the  waveform,  and  k  is  the  nunber  of  first  reconstructed 
points  that  er  will  be  based  upon.  The  waveforms  are  assumed  to 
be  ?56  points  long.  As  seen,  er  is  an  average  square  error  between 
the  filtered  and  reconstructed  waveforms.  When  scanning  is  com¬ 
pleted,  the  window  with  the  smallest  er  is  chosen  as  the  base  for 
reconstruction.  Other  error  criteria  are  possible.  This  criterion 
will  be  used  for  selecting  the  proper  or  best  base  window  to  predict 
the  early  target  response.  Experience  has  proven  it  to  be  effective. 
It  should  be  understood  that  the  reconstructed  tunnel  response  is 
only  true  up  to  the  arrival  time  of  the  actual  tunnel  response. 

The  reconstructed  section  previous  to  this  time  is  just  ficticious, 
but  it  will  be  found  important  for  evaluating  more  precisely  the 
distance  to  the  target.  A  theoretical  study  in  the  next  chapter 
demonstrates  that  the  first  zero  crossing  of  the  reconstructed 
waveform  after  its  deviation  from  the  filtered  corresponds  well 
to  the  target  (tunnel)  -  antenna  distance. 


E.  A  Discussion  on  the  Reconstruct  ion  Interval,  Tr 


Another  important  aspect  to  be  examined  is  the  reconstruction 
interval.  It  must  always  satisfy  the  Nyquist  criterion  for  the 
tunnel  pole  pair  to  be  reconstructed.  It  does  not  play  as  important 
a  role  as  the  pole  extraction  interval.  Tf  should  never  be  too 
small,  for  example  equal  to  T-.  Then,  the  process  may  fail.  It 
will  attempt  to  reconstruct  noise  since  it  is  within  its  frequency 
band.  An  example  is  shown  in  Figure  18.  Here  the  filtered  waveform 
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figure  18.  Example  of  the  reconstruction  process  mhen 
using  a  small  T  . 


in  Figure  17  is  reconstructed  after  bandpass  filtering  through  a 
trapezoidal  filter  with  corners  at  5,100  and  175  Wz.  This  filtering 
removed  any  dc  term  and  high  frequency  components.  The  chosen  Tf 
was  small  ( 4 T g )  as  compared  to  the  tunnel  pole  pair 
[-38.  Meganepers/sec  ♦  j? nx 19.37  Megarad/sec]  Nyquist  interval.  It 
is  seen  that  some  noise  is  carried  along  in  the  reconstructed  section. 

A  reconstruction  interval  close  to  that  required  by  the  Nyquist 
criterion  for  the  tunnel  natural  resonance  is  usually  best.  But 
this  may  require  a  long  base  window.  Such  a  window,  within  which 
the  tunnel  pole  pair  is  described  correctly,  is  usually  not  available. 
For  example,  a  pole  pair  at  20  MHz  has  a  Nyquist  Interval  of  25T_ 
(Tg*?00/?55  nsec).  A  base  window  of  about  50Tg  would  then  be  needed, 
it  was  <i>und,  the  interval  used  by  Prony's  method  for  calculating 
the  natural  resonances  (poles)  works  well  for  the  pole  extraction 
and  reconstruction  processes.  This  interval  was  used  for  recon¬ 
structing  the  filtered  waveform  in  Figure  17.  The  result  is  shown 
in  Figure  19.  Comparing  Figures  18  and  19  we  see  that  use  of  a 
larger  reconstruction  interval  ( N_  =  9 )  reduces  the  reconstruct  ion 
of  noise.  Therefore,  Tr  can  be  effectively  used  for  filtering  out 
undesired  high  frequency  components. 

The  reconstructed  waveforms  are  calculated  by  the  computer 
subroutine  given  in  Appendix  0.  Another  example  of  the  process 
is  shown  in  F igure  20.  The  filtered  waveform  is  the  result  after 
the  extraction  of  the  antenna  pole  pair  [-43.6  Meqanepers/sec 

♦  j2»x69.0  Megarad/ secj  of  the  echo  in  Figure  16(c).  It  was  also 
filtered  by  the  trapezoidal  filter  discussed  previously.  The  re¬ 
constructed  tunnel  pole  pair  was  located  at  [-42.6  Meganepers/sec 

♦  j?wx?1.9  Meqarad/sec}. 

F.  General izat ion  of  the  Reconstruction  Process 


The  reconstruction  process  as  given  in  Equation  (42)  can  only 
reconstruct  one  pole  pair.  The  antenna  pole  pair  and  all  false 
target  resonances  must  first  be  extracted  before  its  application. 
Therefore,  Equation  (42)  can  only  be  used  to  reconstruct  a  target 
response  consisting  of  only  one  natural  resonance  (tunnel).  Most 
targets  are  character ized  by  more  than  one  pole  pair.  For  the  sake 
of  completeness  we  will  Just  derive  the  difference  equation  used 
for  reconstruct ing  target  responses  containing  more  than  one 
resonance. 

Let  us,  for  example  assume  a  target  response  characterized 
by  two  pole  pairs  (complex  conjugate).  After  the  removal  of  the 
antenna  pole  pair  and  other  false  target  resonances  from  its  back- 
scattered  response  then  we  wish  to  reconstruct  the  early  target 
response  for  determining  its  depth.  In  this  case  Equ|t1on  (42)  * 

must  be  modified  to  reconstruct  two  pole  pairs,  (*,,*,)  and  (z,.*?)* 
This  can  be  accomplished  in  a  similar  way  to  that  used  for  deriving 
Equation  (42). 


Figure  20.  Another  example  of  the  reconstruction  process. 


The  7 -transform  of  the  filtered  waveform  (in  its  late  portion) 
can  be  given  as 


Rc(*> 


Y  yJ 

_.±o _ 


where  (Zi.z.)  and  (/•,,/,)  are  the  target  pole  pairs.  After  expansion 
of  the  denominator  ahd  some  manipulation  we  obtain, 

( l-c,z‘  '♦c?**?-c3r3»c4/'4)RcU)*A(2)  (46) 

where  the  constants  c.,  c9,  c-,  and  c,  are  the  same  as  those  in 
Equation  (15).  c  3 

Following  the  same  reasoning  used  for  deriving  Equation  (4?) 
we  can  finally  arrive  at 

rct"VkV4Trl  -  c3rc("Ir*llTB-3rr|- 


Yc<"VkV?V‘circ<"VkVV-cWnVkV  • 


nT  >t  »m  N  T„*4T  ,  k«0,1,...,N  -1 

r-  c  e  e  8  r  *  *  '  '  r 

n«0, 1 , . . . 


where 


c«  .  1  c.  .  ll  rl  .  f*  c.  .  CJ 

co  c4  •  C1  c4  •  c?  c4  •  c3  c4 


and  all  other  parameters  are  as  defined  previously  (see  Equation 
(43)). 

A  generalized  reconstruct  ion  process  analogous  to  Equation 
(18)  (the  generalized  pole  extraction  process)  can  also  be  derived. 
If  a  tarqet  is  character i zed  by  m  single  complex  poles,  then  the 
filtered  waveform  obtained  after  the  extraction  of  the  antenna  or 
other  false  tarqet  resonances  can  be  generally  expressed  as 


(48) 
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where 


r 

r- 


L  Wv 

r 


m 


m 


Also,  note  again  that  in  the  above  equation  mr  does  not  represent 
the  target's  complex  conjugate  pole  pairs  but  just  its  single  poles 
which  are  twice  in  number  than  the  complex  conjugate  pole  pairs. 

As  a  check  one  can  easily  derive  Equations  (42)  and  (47)  from 
the  generalized  Equation  (49).  It  should  be  understood  that  con¬ 
current  reconstruction  of  a  large  number  of  poles  is  Impractical 
since  it  requires  a  large  base  window  for  the  process  to  begin. 

Also,  note  that  any  dc  term  must  be  removed  before  applying  Equation 
(49). 

G.  Assets  of  the  Reconstruction  Process 

The  reconstruct  ion  process  is  an  important  part  of  the  tunnel 
identification.  It  accomplishes  three  tasks: 

1.  The  clutter  problem  is  solved  as  long  as  it  appears  in 
the  early  portion  of  the  echo  (it  usually  is).  Exact  knowledge 
of  the  tunnel  response  is  not  needed.  Only  a  small  window  with 
fair  accuracy  is  required  so  the  true  tunnel  response  can  be  isolated. 
The  reconstruction  process  will  predict  the  rest  of  the  waveform 
based  on  this  window. 

?.  Host  of  the  times  the  later  portion  of  the  received  echo 
is  not  noise  or  clutter  free.  Based  on  the  error  criterion  in 
Equation  (44)  the  reconstruction  process  can  detect  the  window  which 
best  describes  the  tunnel  pole  pair.  Consequent ly,  a  more  accurate 
early  time  response  can  be  derived.  Equation  (44'  is  minimized 
when  uniformity  occurs  in  the  base  window.  If  noise  or  clutter 
is  present  in  the  base  window  under  testing  the  predictability  of 
the  reconstruction  process  Is  disrupted,  this  causes  an  error  in¬ 
crease.  Note,  also,  that  use  of  the  error  criterion  satisfies 
Equation  (43). 

3.  In  the  next  chapter  it  will  be  demonstrated  that  the  re¬ 
constructed  waveform  can  play  a  very  important  role  in  the  identi¬ 
fication  process.  A  comparison  of  the  filtered  and  reconstructed 
waveforms  can  give  us  a  very  qood  estimate  of  the  tunnel’s  response 
arrival  time.  Thus,  its  distance  from  the  antenna  can  be  estimated. 

In  turn,  mapping  can  present  a  complete  picture  of  the  tunnel's 
depth,  structure  and  relative  position. 
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!n  the  next  chapter  a  simple  experimental  model  of  the  Radar- 
Tunnel  Structure  is  examined.  The  study  uses  the  pole  extraction 
and  reconstruction  processes  in  order  to  extract  information  from 
the  theoretical  echoes  and  determine  the  limitations  involved. 

In  Chjpter  V  the  results  of  the  model  are  used  and  a  study  on  measured 
echoes  is  pursued. 


r 
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CHAPTER  IV 

STUDY  Of  A  S1MPLF  TRANSMISSION  LINE  MODEL  OF  AN 
UNDERGROUND  RADAR -TUNNEL  STRUCTURE 


A.  Object  ives 

The  (jo,ih  of  this  chapter  are  the  following: 

1.  To  present  a  physical  and  mathematical  study  of  a  trans¬ 
mission  1  me  model  which  approximates  the  actual  tunnel  scattering 
behavior. 

?.  To  relate  the  model  to  the  actual  radar-tunnel  structure 
and  derive  methods  to  extract  information  from  the  theoretical 
echoes.  These  methods  can  then  be  used  when  dealing  with  measured 
echoes. 


3.  To  demonstrate  the  role  and  significance  of  the  pole 
extraction  and  recons troct ion  processes  to  be  used  in  the  identi 
f icat ion  process. 

B.  The  Model 


We  are  attempting  to  construct  an  experimental  model  of  the 
Radar-Tunnel  structure  to  study  and  extract  information  from  the 
received  backseat tered  echoes.  This  will  help  us  to  analyze  the 
measured  responses. 

The  input  to  the  antenna  must,  of  course,  have  a  wide  spectrum 
to  ensure  excitation  of  the  tunnel's  natural  resonances.  It  is 
unfortunate  though,  that  this  is  not  easily  achieved  since  the  dipole 
antenna  is  in  itself  a  resonant  structure.  If  the  input  has  a  wide 
spectrum,  say  a  gaussian  pulse,  the  antenna  will  act  as  a  filter. 

For  excitation  the  tunnel’s  natural  resonance  must  occur  close  to 
the  antenna’s  resonance. 

Since  the  target  and  the  antenna  are  both  resonant  structures, 
they  can  be  modeled  as  low-Q  ana’og  filters.  Such  can  be  a  simple 
RLC  circuit  shown  in  Figure  ?1.  The  impedance,  7,  of  the  network 
is  given  by 


s?/C 


s'cs/Rf  ♦  1/LC 


(SO) 


5? 


The  poles  of  Z  are  located  at 


s.s  •  *  yip  !  J ,  ft  *  (r^) 


(51) 


The  theoretical  model  can  then  be  constructed  of  two  RtC  net 
works  separated  by  a  lossy  transmission  line  of  characteristic 
impedance  ZQ(.>).  The  transmission  line  is  introduced  to  simulate 
ground  medium  to  some  degree. 

The  input  to  the  system  will  be  a  gaussian  pulse  generator 
w'd“  spectrum'  of  internal  resistance  R,.  as  shown  in  Figure  22. 


GAUSSIAN  PULSE 
GENERATOR 


Figure  22.  Transmission  line  nodtl. 
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’••rviances  ZA  and  Zy  of  the  model  represent  the  antenna  and  target 
(tunnel),  respect  i  vely.  They  are  of  the  form  of  Figure  ??.  ZA 
will  consist  of  resistance,  Rft,  capacitance,  C^,  and  inductance, 
lA.  Similar  definitions  will  stand  for  Ry,  Cy  and  Ly  of  Zy. 

Before  proceeding  to  the  mathematical  analysis  of  the  model, 
let  us  first  give  it  a  physical  interpretat ion.  For  convenience 

the  me.iel  is  redrawn  as  in  F  igure  ?3. 


Figure  ?3.  Transmission  line  model  indicating 
wave  propagation. 


The  output  V  (s-j^i)  in  Figure  ?3  will  represent  the  response 
of  the  system.  A  gauss lan  pulse  is  used  as  a  source  at  t*0. 

2  2. 

(A  gauss i an  pulse  is^defingd  as  g(t)  •  e'**’*’*  '  '*P,  where  t  is  its 
time  width.)  At  t»0*.  .a  voltage  appears  at  R  and  also  across  ZA,(a). 
This  is  shown  in  Figure  ?4  where  the  forced  response  and  natural 
response  of  the  circuit  are  illustrated.  The  transmission  line 
is  also  excited  and  a  waveform  similar  to  Figure  ?4  propagates  toward 
Z..  As  it  reaches  the  end  of  the  transmission  line,  Zy  is  excited 
and  a  reflected  wave,  (c),  is  now  launched  back  toward  the  generator. 
It  contains  resonances  associated  with  both  7.  and  Zy.  This  reflected 
waveform  reaches  Z.,  and  creates  an  output  voltage  (v  (s)),  (e). 

If  the  parallel  combination  of  R  and  7.  Is  not  matched  to  Z  ,  then 
a  portion  of  it  is  also  reflected  back  toward  Zy,  and  the  process 
continues.  A  lonq  time  waveform  is  shown  in  Figure  ?5.  Here,  we 
will  study  the  general  case  in  which  the  parallel  combination  of 
R  and  Z.  is  not  matched  to  Z  ,  and  discuss  its  limitations  as  the 
transmission  line  becomes  short. 
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Figure  ?4.  Response  of  an  RlC  resonant  network  to  a  Gaussian  Input  pulse 


Our  goal  is  to  isolate  the  resonances  associated  with  Zy  and 
to  evaluate  the  length  of  the  transmission  line.  The  procedures 
introduced  in  the  first  three  chapters  are  applied  to  isolate  the 
desired  parameters. 


Since  we  already  know  them  precisely,  we  can  then  better 
understand  the  limitation  of  our  techniques  when  applied  to  unknown 
systems,  such  as  a  tunnel.  Observe  that  the  data  will  be  corrupted 
by  the  presence  of  whatever  parameters  are  used.  One  set  of 
parameters  is  studied  here  in  detail,  but  results  for  various  other 
parameters  are  given  in  Appendix  E. 

C.  Mathematical  Analysis  of  the  Model 


A  quantitative  description  of  the  transmission  line  model 
will  now  follow. 

The  method  of  analysis  pursued  in  this  study  involves  trans¬ 
forming  the  gaussian  pulse  (generator  output)  to  the  frequency  domain, 
evaluating  V  (jj  and  then  inverse  transforming  it  to  the  time  domain, 

(vQ(t)).  0 

Two  approaches  will  be  used  and  compared  for  calculating  the 
system's  response  of  Figure  2 3.  The  first  follows  the  actual  physical 
multiple  reflection  process  just  described  and  transforms  the  desired 
voltage  to  time  domain  after  each  reflection.  The  second  obtains 
the  avnplote  solution  in  the  frequency  domain  and  transforms  the 
total  voltage  to  time  domain.  The  approaches  are  equivalent  and 
both  are  used  here  to  ensure  that  no  computational  errors  are 
generated. 


Defining  the  input  as  a  gaussian  pulse  of  spectrum  G(s*Jw), 
the  propagation  constant  as  y*o*Jf,  and  reflection  and  transmission 
coef f ic ieots  as 


V*> 


‘Ao 


(52a) 


Z.  -Zrt 
i  \  AS  O 

pA{s)  *  77TT 

As  0 


?  7. 

,  .  As 
’A'  ’  ’  ^0 


(52b) 


(52c) 


vh-*"-  ■ 


the  voltdqe  at  the  indicated  points  in  Figure  ?3  can  be  expressed 


as  follows: 

Va(s)  *  G(s)  Ty(s)  (53a) 
V„(s)  «  G(s)  tv(s)  (53b) 
Vc(s)  *  G(s)  ty(s)  Pj(s)  e'Yd  (53c) 
vd(s)  •  G(s)  Tv(s)  ^(s)  e’?Yd  (53d) 
Ve(s)  •  G(s)  ty(s)  pj(s)  tA(s)  e‘2rd  (53e) 


V  (s)  is  the  first  echo  seen  across  Z A-  Consequently,  the  first 
••cho  across  can  be  expressed  as 

VfSt(s)  •  G(s)-Ve(s)  *  -  Ve(s).  (54) 

It  appears  after  a  tine  delay  of 


or 


t 


?d 

TnTtl 


t  *  ?d/cii  «  6.667  d/FJi IT  (nsec)  (55) 

for  lossless  media,  d  is  in  meters,  and  c_  and  u_  are  the  relative 
permittivity  and  permeability  of  the  medium, respect ively. 
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The  second  echo  across  R  will  be  delayed  by  13.333  cV' c rw*r 
nsec.  According  to  the  above  analysis  it  is  given  by 

V„nd(s)  *  F-'v(s)  ^(S)  nA{s)  e*4Yd]  G(j,j)  *  {56) 

The  continuous  time  voltage  seen  across  R^  is  the  sum  of  all 
the  echoes,  J.e. , 


*„<*>  •  v>'5)  *  Vs’  *  V',(5)  ‘  - 


(57) 


In  obtaining  the  frequency  response  of  V  (jw)  we  substitute 


\(s) 


ta(s) 


in  Equations 

(52). 
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(58b) 


(58c) 
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The  poles  of  t  (s),  o.(s)  anJ  T.(s)  are  located  at 

V  A  t\ 


(58d) 
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The  imaginary  parts  of  Equations  (59)  and  (60)  correspond 
to  the  complex  resonances  of  the  coefficients  of  Equations  (5?). 

Furthermore,  if 


then  the  complex  resonances  of  the  coefficients  will  be  approximately 
equal  to  the  natural  resonances  of  Z.  and  Z,,  as  defined  in  Equation 
(51).  If  Equation  (61)  is  not  satisfied  then  the  resonances  of 
the  coefficients  will  deviate  from  the  actual  ones  of  the  target 
or  the  antenna,  and  the  target  structure  can  not  be  determined 
correctly  according  to  the  relationships  given  later.  In  all  Our 
studies  we  will  always  attempt  to  satisfy  Equation  (61). 

According  to  the  above  discussion  the  first  echo,  V  _,(s), 

*  lsl  *. 

has  double  pole  pair  at  (s*-.sa>)  single  pole  pair  at  (Sj.,,Sy) 

Consequently,  V  .(s)  is  characterized  with  triple  pole  pair  at 

•  7™  * 

(s.,,sA,)  and  double  at  (s~,,s*,).  The  poles  of  the  subsequent 
ecnoesxan  be  determined  in  a  Similar  way. 

The  above  analysis  of  the  transmission  line  model  can  serve 
the  purpose  of  a  good  understanding  of  the  response,  but  it  becomes 
difficult  to  implement  for  a  computer  analysis. 
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A  more  concise  method  of  looking  at  the  transmission  line 
is  to  represent  an  equivalent  impedance  of  Zy  at  the  input  terminals 
of  the  line.  According  to  the  wave  equations  governing  the  trans¬ 
mission  line,  the  equivalent  impedance  of  Zy  (Figure  22)  at  the 
input  terminals  is  given  (?7)  by 


Z0*Zytanh(yd) 

^Teq  Zo  7y*Z0tanh( ^dT  * 

Representing  the  parallel  combination  of  ZTeq  and  Zft  as 

J ,  . 

<"  ^te.,  - 


then  VQ(s)  can  be  expressed  as 


\  s  in 

VQ(s)  •  W - JTJ-  *-•  -  B 

0  Zin*Rs  VZin 


The  above  expression  is  equivalent  to  that  given  in  Equation  (57). 
0.  Testing  of  the  Model 

Equation  (64)  was  used  for  implementation  of  a  computer  model 
for  the  transmission  line,  given  in  Appendix  E.  The  characteristic 
impedance  of  the  line  in  the  computer  model  was  defined  as 


-  . 

o  >1 


and  the  propagation  constant  as 


r*juivJruc 


The  parameters  r  and  o  represent  the  series  and  shunt  losses  per 
unit  lenqth  of  the  line,  while  u  u  and  crc  correspond  to  the 
inductance  and  capacitance,  repsecfi vely,  per  unit  length  of  the 
1  ine. 
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It  was  necessary  to  Introduce  some  loss  in  the  line  (complex  y) 
in  order  for  v  (t)  (the  inverse  Fourier  transform  of  V  (jw))  to 
be  completely  diminished  within  the  chosen  time  window,  T  .  A  finite 
time  window  is  a  limitation  of  the  Fast  Fourier  Tr ansf ortn*used  for 
time  inversion  of  the  819?-point,  complex,  discrete  array,  VQ( J^.) . 
This  array  was  produced  by  an  incremental  frequency  interval  equal 
to 


\>j  *  (rad).  (67) 

w 

The  above  matters  should  be  carefully  arranged  for  the  successfulness 
Of  the  computer  model. 

Furthermore,  for  simplicity,  2  was  kept  real  by  defining 

r  as 


r 


vv 

cOcr 


(68) 


The  output  voltage  shown  in  Figure  ?$  was  calculated  by  the 
computer  model  in  Appendix  E.  The  input  was  a  6  nsec  gaussian  pulse. 
The  pole  pair  of  was  placed  at 

^A.s^  •  -  ?0x  I0fi  nepers/sec  ♦  J?m50x  10^  rad/ sec 


and  that  of  Zj  at 

Sj.Sj  *  -  10x10^  nepers/sec  ♦  J?*x4QxlO^  rad/sec  . 


Other  parameters  as  indicated  in  Figure  ?5  are 


ra  *  ?son 
rt  •  soon 
Z0  •  J77Q 
Rs  «  3770 
d  *  30  meters 


8.84S6X10*1'  Farads/meter 
1.766x10*^  Henrys/meter 
5xl0"fi  tVmeter 
1. 

1. 
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Not*  *  i it  1  was  chosen  long  enough  so  that  there  would  be  no  overlap 
of  the  echoes.  The  pole  pair  of  Z.  was  selected  with  respect  to 
that  of  Z  so  ihat  the  maximum  response  of  Z.  alone  in  the  frequency 
hvnain  would  occur  at  the  frequency  where  th£  response  from  7^ 

(alone)  would  decay  to  its  half  power  value. 

As  ment ioned  in  Chapter  I,  we  are  interested  in  determining 
the  structure  of  the  tunnel  (target)  by  examining  the  received  voltage. 
The  natural  resonance  of  the  tunnel  (imaginary  of  sT)  will  determine 
its  structure.  The  peak  response  from  the  tunnel  occurs  when  the 
reflection  from  the  front  and  back  interfaces  are  in  phase.  To 
a  f irst  order  approx imat ion,  neglecting  caustics,  this  occurs  for 

?BM  -  «  ,  (69a) 

where  H  is  the  tunnel  height.  Assuming  the  velocity  of  the  wave 
in  the  tunnel  to  be  that  of  free  space.  Equation  (69a)  gives 

H  •  (69b) 

H 

where  f  is  a  comparable  first  order  resonance  corresponding  to 
the  height  of  the  tunnel. 


Its  depth,  d.  can  be  determined  by  the  arrival  time,  T  ,  of  the 
first  echo  as  follows: 

T  (nsec) 

d  - - - - -  (meters).  (70) 

6.667  •  t 


From  the  above  discussion  it  is  instructive  to  use  the  first 
echo  for  analysis.  Since  this  is  the  simplest  response  containing 
information  about  7,.  The  time  expanded  first  echo  of  Figure  ?5 
is  shown  in  Figure  c6(»).  For  comparison,  the  first  echo  as  calculated 
by  Equation  (54)  is  also  shown  in  Figure  ?6(b).  Note  that  it  is 
identical  to  that  in  Figure  ?6 (a).  The  poles  of  this  echo  were 
found  by  Prony’s  method  to  be  located  at 


tl  1<5A  001  t  (ll 


p,.p, 


18.5705  Me<ja»epers/sec  ♦  >  7 . OS08  M-egar ad/sec 


P?.P?  •  -  47.09545  Meganepers/sec  1  j?»x60. 67961  Megarad/sec 
P3.P3  *  73.73658  Meqanepers/sec  *  ,P^>?9. 83078  Megarad/sec 

with  respective  residues  of 

Rp  .R*  -  (2.866151  ♦  J7. 846461 )xlO'? 

Rp  .R*  -  (-7.354673  ♦  j3.4901?9)x10‘? 

Rn  ,R*  •  (-.4535518  ♦  J.07051006)xl0‘?  . 

pt 

( The  residues  when  calculated  by  Prony's  method,  do  not  have  any 
absolute,  but  only  relative  roeaninq,  since  they  depend  on  the 
Starting  point  of  the  window  used  by  Prony's  method  for  finding 
the  poles.)  If  Equations  (5?)  are  used  for  calculating  the  poles 
°f  V  (s),  we  obtain 
l'1 

Pj.P'j  *  -  46.575  Meqanepers/sec  ♦  .17**49.551  Meqarad/sec 

Pp.p'i  •  -  46.575  Meqanepers/sec  *  j?w*49.551  Megarad/sec 

pj.p* 3  1  -  73.763  Meqanepers/sec  ♦  J?«x39.861  Megarad/sec 

These  results  agree  very  well  with  the  ones  generated  by  Prony's 
method.  Note  that  the  tunnel  height  corresponding  to  a  39.83  MH* 
resonance  is  found  to  be 


H  •  - c- - -  .  1.883  m. 

4x39.83x10° 


Tlic  irrival  time  in  this  case  is  easily  detected  to  be 
T.  «  6.667(30)  •  POO  nsec 

a 

since  the  output  data  were  produced  from  a  theoretical  model. 

E .  Use  of  the  Reconstruct  ion  Process  to  Estimate 
!he  Arrival  Time  of  The  Target  Response 

It  was  noted  previously  that  in  general  when  actual  echoes 
are  encountered  the  early  portion  of  the  echo  is  dominated  by  clutter. 
Thus,  we  will  attempt  to  describe  a  more  clutter  immune  method  for 
calculating  the  arrival  time  of  the  target  response. 

Employing  the  pole  extraction  process,  the  antenna  double 
pole  was  extracted  from  the  response  of  Figure  ?6.  The  result  is 
shown  in  Figure  ?7.  Note  the  peaking  at  the  beginning  c*  t‘  re¬ 
sultant  waveform.  This  is  because  the  forced  response  ri  ’  be 
characterized  with  a  set  of  poles.  The  rest  of  the  Jilt  >veform 

contains  only  the  pole  pair  of  Zy  (corrected),  (Pj.Pj)- 

After  removing  the  dc  component  of  the  target  respc  r  t^u'e 
?7),  we  can  then  reconstruct  its  early  portion  based  on  it  late 
response,  as  shown  in  Figure  ?8.  Concentrat ing  in  the  inter. al 
of  19-37  nsec  we  observe  that  in  this  time  span  the  reconstructed 
waveform  deviates  drastically  from  the  course  of  the  filtered  wave¬ 
form.  The  first  zero  crossing  of  the  reconstructed  portion  after 
its  deviation  from  the  filtered  waveform  corresponds  to  the  tail 
end  of  the  forced  response  (see  Figure  ?6).  This  phenomenon  was 
tested  for  many  such  responses  in  order  to  examine  its  validity 
(see  Appendix  E).  It  was  found  to  be  the  same  in  every  case  that 
was  studied.  From  Figure  ?8  the  target's  depth  can  be  calculated 
as 


>1 


6 . 66  7  ^trvr 


?00 

575T7 


30  meters. 


(71) 


where  t  is  the  gauss ian  pulse  width  in  nsec  and  T  is  the  time  in  nsec 
correspBnding  to  the  first  zero  crossing  of  the  reconstructed  wave¬ 
form  after  its  deviation.  (Note  that  the  origin  of  Figure  ?8  cor¬ 
responds  to  180  nsec,  see  Figure  ?6.) 

The  previously  studied  response  was  obtained  with  a  long  trans¬ 
mission  line  (3Qn}.  The  echoes  then  did  not  overlap.  We  can  now 
shorten  the  line  in  order  to  determine  if  it  is  still  possible  to 
detect  the  target.  By  shortening  the  line,  there  will  be  overlap 
of  the  echoes  and  probably  the  antenna  response,  V  ( s )  (Equation 
(53a)).  We  will  neglect  the  antenna  response  since  for  the  parameters 
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traction  of  antenna  double  pole  pair  from  first  echo. 


-6  NSEC  TIME  REFERENCE:  0*180  NSEC 
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chosen  it  Mould  be  negligible  by  the  arrival  of  the  first  echo. 

If  it  did  not,  then  an  increased  residue  of  the  antenna  pole  pair 
Mould  have  to  be  encountered  in  the  first  echo. 

Using  the  same  parameters  as  those  for  dOO«n.  on  page  65, 
the  f  irst  echo  response  for  depths  ?0,  18,  15  and  10  meters  are 
shown  in  Figures  29,  30,  31  and  3?  (dotted  lines),  respectively. 

Note  the  Mindow  in  each  Maveform  containing  the  first  echo  only. 

This  Mindow  was  used  by  Prony's  method  for  calculating  the  poles 
of  the  first  echo  (sampling  interval  is  Tg*6000/8191  nsec).  We 
could  not  use  the  total  window  shown  in  tRe  figures  since  Prony's 
me t hod  assumes  that  the  response  associated  with  all  the  poles  in 
the  region  being  sampled  has  been  introduced  pior  to  the  original 
sample. 

For  example,  let  us  randomly  select  a  time  window  out  of  the 
total  span  of  the  waveform  and  use  Prony's  method  to  find  its  poles. 

In  this  time  window  there  could  be  a  superposition  of  two  or  more 
echoes.  It  was  found  that  Prony's  method  will  probably  find  all 
of  the  comple*  resonances  of  the  poles  but  the  calculated  real  parts 
will  be  far  off  from  the  ones  existing  in  the  waveform  if  a  resonance 
is  initiated  within  this  time  window.  In  fact  Prony's  method  re¬ 
quires  that  no  forced  response  be  present  in  the  data  if  the  real 
parts  of  the  poles  are  to  be  calculated  accurately. 

The  poles  obtained  by  Prony's  method  for  the  echoes  in  Figures 
29,  30  and  31  are  shown  in  Table  I.  These  are  much  the  same  with 
those  obtained  for  d»30m.,  as  expected.  In  order  to  determine  the 
depth  of  the  target  we  first  extract  the  antenna  double  pole  pair. 

The  result  is  the  solid  line  on  the  same  f iqures  (Figures  29,  30, 
and  31).  After  centering,  the  reconstructed  early  portion  of  the 
target  response  is  shown  in  Figures  33,  34  and  35  for  depths  at 
20,  18  and  15  meters,  respectively.  Detecting  the  first  zero  crossing 
of  the  reconstructed  target  responses  after  their  deviation  from 
the  filtered  waveform  (solid  line)  we  obtain  the  following  results: 


Tq  »  120  ♦  20. 12  •  140.2  nsec  ;  d«20m. 

T^  *  108  ♦  17.8  •  125.8  nsec  ;  d*l8m. 

Tq  •  90  *  16.5  »  106.5  nsec  ;  d«l5m. 

Using  the  above  results  in  Equation  (71)  the  depths  are  calculated 
to  be  20.13,  17.97,  15.07  meters  as  compared  to  the  actual  depths 
of  20,  18  and  15  meters,  respect i vely. 
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Figure  30.  Extraction  of  antenna  double  pole  pair  from  first  echo. 
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antenna  double  pole  pair  from  first  echo  (d=10m.) 


fable  1 

PRONY  RESULTS;  MIN.  SQUARE  ERROR  CASE 


‘WVFRM 

DEPTH 

IBS1 

— 

ITS? 

P0LES3(k106) 

RESIDUES(xl0'?) 

REAL 

IMAG 

REAL 

IMAG 

Fig.  ?6 

30  m. 

6lJ 

SI 

-43.5?05 

-42.09S4 

-23.23658 

♦47.95888 

♦50.62961 

♦39.83028 

2.866151 

-2.354673 

-  .4535518 

♦2.846461 

♦3.490129 

♦  .07051006 

Fig.  ?9 

20  m. 

5Tb 

4S 

-S3. 0749 

-40.19353 

-25.16282 

♦47.19403 

♦50.76860 

♦40.15106 

.1067381 

.8896122 

-  .4398645 

♦2.831341 

♦2.994219 

♦  .3466342 

Fig.  30 

19  m. 

5Tb 

44 

-43.58974 

-44.69857 

-25.06127 

♦47.65909 

♦51.03842 

♦39.62526 

.842728 

.0540642 

-  .2190129 

♦3.228905 

♦3.139947 

♦  .408473 

Fig.  31 

IS  m. 

4Tv 

40 

-35.51768 

-54.8633 

-21.13081 

♦43.77227 

♦51.44039 

♦39.51648 

2.961494 

-2.166434 

-  .301001 

♦  .2245375 

♦  .5655805 

♦  .212935 

Fig.  3? 

10  m. 

First  Echo  Time  Window  not  Wide 

Enough  for  Processing 

*IBS  *  Interval  between  samples  used  by  Prony's  method  for 
calculating  the  poles. 

ITS  -  Starting  point  of  the  window  used  by  Prony's  method 
for  calculating  the  poles. 

The  real  and  imaginary  parts  of  the  poles  are  given  in 
Nepers/sec  and  Hz,  respectively.  The  first  two  pole  pairs 
correspond  to  the  antenna  double  pole  pair  and  the  third 
to  the  tunnel. 

\  *  .73251  nsec  (6000/8191) 
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As  seen,  the  reconstruct  ion  process  has  the  capability  of 
producing  very  accurate  results.  This  is  provided  the  filtered 
waveform  does  contain  a  time  window  where  the  target  response  is 
actually  existing  by  itself. 

Let  us  now  turn  our  attention  to  Figure  3?  (dotted  line). 

In  this  case  the  line’s  length  was  at  10  meters.  The  multiple  echoes 
are  66.66?  nsec  apart  and  they  are  overlaped.  Also,  the  arrival 
time  of  the  first  echo  (66.667  nsec)  Is  short  enough,  so  the  antenna 
(.'»)  rtspcrsc  interferes  with  the  first  echo's  response.  Therefore, 
the  time  window  containing  the  first  echo  only,  as  shown,  is  not 
wide  enough  in  order  for  Prony's  method  to  determine  the  poles  of 
the  first  echo.  We  could  have  superseded  this  problem  if  our 
samp  1 ing  interval,  Tg,  was  smaller.  Thus,  more  processing  points 
would  be  available  in  the  respective  time  window.  These  points 
could  then  be  enough  for  generating  the  equations  required  by  Prony's 
method  to  calculate  the  poles.  But  this  would  require  larger  array 
for  v0(j«/.  Note  that  Tg  is  given  by 

T 

8  *  ITT 

where  K  is  the  n^nber  of  points  describing  the  complex  array  of 
V  (j  ).  The  computer  used,  could  only  process  a  maximum  of  <*819? 
iRmplei  points,  which  were  used  for  all  the  data  presented  here  and 
in  Appendix  F. 

This  seems  to  be  the  limit  of  our  processing  capability. 

As  the  line  becomes  shorter  more  of  the  multiple  echoes  overlap 
and  the  waveforms  become  more  difficult  to  process.  Its  main  compli¬ 
cation  is  again  that  poles  are  introduced  after  the  starting  point 
of  the  minimum  time  window  required  for  Prony's  method  to  calculate 
all  the  poles  present  in  that  window.  If  we  are  to  accurately  process 
a  response,  the  chosen  time  window  must  be  coherent  and  continuous. 

Since  the  pole  of  the  first  echo  are  known  *r am  previous 
processing  we  can  go  one  step  further  with  the  response  at  d«10m. 

Usinq  the  antenna  double  pole  pair  as  obtained  for  d*30m. ,  it  is 
extracted  from  the  echo  of  Figure  3?  (dotted  line).  The  result 
is  the  solid  line  in  the  same  figure.  Note  the  peaking  of  the  forced 
response.  They  occur  at  the  beginning  of  each  echo,  where  the  forced 
response  is  confined.  These  peaking  show  that  three  echoes  exist 
in  the  total  time  interval  of  the  waveform. 

The  target  response  alone  is  Indicated  in  the  first  echo  for 
a  time  length  corresponding  to  a  cycle.  Note,  also,  that  the  target 
response  is  not  isolated  from  the  other  echoes,  yet,  since  they 
consist  of  more  poles.  This  was  developed  earlier.  Using  the  target 
pole  pair  for  d»30m.,  this  target  response  alone  was  reconstructed 
as  shown  in  Figure  36.  The  first  zero  crossing  of  the  reconstructed 
waveform  after  its  deviation,  is  detected  to  be 
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Reconstruction  of  target  response  alone  (d*30m.). 
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7^  n\ec. 

Therefore, 


d  *  57557  '  ,0-05  m- 

Indeed,  this  is  the  actual  length  of  the  transmission  line. 

The  above  observation  indicates  that  if  a  portion  of  the  target 
response  exists  in  a  small  time  window  then  the  reconstruction  process 
can  still  be  used  for  detecting  the  target's  distance  from  the  antenna 
(or  constructing  a  map). 

F .  Cone  1  us  ions 

It  is  concluded  that  Our  main  problem  with  processing  the 
echoes  is  the  detenu inat ion  of  the  poles  by  Prony's  method.  Once 
the  poles  of  the  waveform  are  known,  within  a  reasonable  accuracy, 
the  pole  extraction  process  can  always  be  used  for  isolating  the 
target's  pole  pair.  In  turn,  the  early  target  response  is  recon¬ 
structed.  The  first  zero  crossing  of  the  reconstructed  waveform 
after  its  deviation  from  the  filtered,  can  be  used  in  Equatton  (71). 
The  target- antenna  distance  is  subseguently  estimated. 

When  the  tarqet’s  depth  is  small  then  the  multiple  reflections 
of  its  response  could  overlap.  If  the  time  window  where  the  first 
echo  appears  without  overlaps  is  not  large,  (greater  than  ?N  or 
T’i  time  points,  where  N  is  the  ntxnber  of  sinqie  poles  in  tne 
window)  then  Prony' s  method  will  fail  to  calculate  the  correct  poles. 
This  can  be  alleviated  by  matching  (see  Equation  (5?))  with 

the  line's  characteristic  impedance,  l  .  Thus,  all  echoes,  except 
the  first,  are  suppressed.  In  our  field  study  at  Gold  Hill,  Colorado, 
discussed  in  the  next  chapter,  the  antenna  structure  was  closely 
matched  to  the  qround  impedance  to  avoid  multiple  tunnel  reflections. 
Of  course,  ground  impedance  is  a  dynamic  parameter.  A  study  could 
tie  attempted  here  to  determine  the  tolerance  of  the  antenna  structure 
and  ground  impedance  mismatch  on  our  results. 

Furthermore,  we  should  not  neglect  that  Equation  (61)  should 
be  reasonably  satisfied  at  all  times.  Otherwise,  erroneous  results 
could  be  obtained  for  the  tunnel  or  target  structure  (see  Appendix 

n. 
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CHAPTER  V 

APPUCAUON  Of  THE  POLE  EXTRACTION  ANO 
RECONSTRUCTION  PROCESSES  ON  MEASURED 
TUNNEL  ECHOES 


A.  Object  Ives 

The  objectives  of  this  chapter  are  the  following: 

1.  To  apply  the  pole  extraction  and  reconstruction  processes 
on  measured  tunnel  echoes  for  detection  and  identif ication. 

?.  To  evaluate  our  results  and  compare  them  with  results 
from  other  processing  attempts. 

B.  Echo  Recording  and  Processing 

During  the  summer  of  1978  a  group  from  the  E lectroSc ience 
Laboratory  recorded  a  set  of  echoes  on  a  tunnel  site  at  Sold  Hill, 
Colorado.  The  position  of  the  recorded  waveforms  with  respect  to 
the  tunnel  is  indicated  in  Figure  37.  As  shown,  two  traverses  were 
made  over  the  rectangular  tunnel  [??]  .  At  position  RIO  the  tunnel 
was  about  20  ft.  (6.1  m.)deep,  and  of  si/e  4'x9'.  These  data  would 
correspond  to  a  tunnel  response  arrival  time  cf  about  100  nsec  (the 
ground  relative  permittivity  was  accurately  measured  to  be  c  "6) 
and  a  resonance  at  about  ?5  Wiz. 

A  6  nsec  gaussian  pulse  was  used  for  excitation  with  an 
8-foot  crossed  dipole  antenna,  referred  to  as  LBANT  (Long  Box  Antenna) 
(14).  Some  of  the  recorded  echoes  are  shown  in  Figure  16.  The 
later  portion  of  all  the  waveforms  was  used  by  Prony's  method  for 
calculating  their  poles.  Results  are  given  in  Tables  2  and  3. 

Two  complex  conjugate  pole  pairs  are  indicated.  The  one  at  40-50 
*</  corresponds  to  the  antenna  pole  pair,  and  the  other  at  ?0-30 
Wz  to  the  tunnel  pole  pair.  The  real  part  of  the  calculated  poles 
varied  over  a  wide  range.  This  is  due  to  the  nonuniformity  of  the 
ground,  and  also  to  the  sensitivity  of  the  real  part  in  the  calcu¬ 
lations  performed  by  Prony's  method. 

The  antenna  pole  pair  given  in  the  tables  was  extracted  from 
each  particular  echo  by  employing  Equation  (16)  and  was  subsequently 
•rrected  for  distortions  according  to  Equation  (34).  The  early 
portion  of  the  tunnel  response  for  each  waveform  was  then  reconstructed 
according  to  Equation  (4?),  with  the  tunnel  pole  pair  as  given  in 
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Figure  37.  Measurement  coordinates  of 
Gold  HUl  Map  2. 


the  tables.  The  reconstructed  waveforms  at  positions  008  and  RIO 
are  shown  in  Figure  19  and  20,  respectively.  It  should  be  mentioned 
that  the  intervals  used  during  the  pole  extraction  and  reconstruction 
processes  were  the  same  as  those  used  by  Prony's  method  for  calculating 
the  waveform  poles.  These  intervals  are  given  in  Tables  2  and  3 
under  the  column  of  I8S  (Interval  Between  Samples). 
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JNY  RESULTS; 


Table  2 
MIN.  SQUARE 


El _ _ 


_ POLES _ RESIDUES _ 

WVFRM 

POS. 

ITS* 

(Nepers* l06/sec,MH;) 3  (x5x!0'3V) 

REAL 

I  MAG 

real 

I  MAG 

MOO 

9Tb 

107 

-52.57866 

-24.8865? 

♦67.33230 

♦25.85206 

.4162106 

. 1878496 

NO? 

9T8 

10? 

-10.79741 

-32.04296 

♦51.30669 

♦24.91175 

-  .2334531 

.0342021 

♦  .1110562 

♦  .2627557 

a 

9Tb 

105 

-17.33083 

-24.33710 

r*47. 07857 
♦20.47811 

.397899? 

-  .702791 

♦  .7173292 

♦  .3943723 

P06 

8Tb 

11? 

-32.70187 

-30.30265 

♦4 1 .9480? 

t»8. 21603 

-  .1942287 

.6458865 

♦1.310707 

♦  .08358793 

Q08 

9Tb 

101 

-26.92599 

-38.10888 

♦45.48387 

♦19.3685? 

.04361706 

-1.484128 

♦  .5112753 

♦  .7730645 

RIO 

7Tb 

96 

-43.5940? 

-  4 . 259986 

♦69.02293 

♦21.8921 

-  .000748475 

-  .5154413 

*  .178899 

♦  .3856364 

SI? 

_ 

98 

-57.71117 

-18.19595 

♦46. 7993? 

♦25. 1 797? 

-  .7286549 

.06751176 

♦  .6220561 

♦  .09678334 

1 1 8S  »  Interval  between  samples  used  by  Prony's  method 
for  calculating  the  poles  (Tg»?00/?55  nsec). 

ITS  *  Starting  point  of  the  window  used  by  Prony's 
method  for  calculating  the  poles. 

^  The  first  pole  pair  corrrsrcnds  to  the  antenna 
and  the  second  to  the  tunnel. 
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Table  3  '  > 

PRONY  RESULTS;  MIN.  SQUARE  ERROR  CASE 


WVFRM 

POS. 

IBS7 

Tt? 

( Nepers* 106/sec,MH/)" 

1  (*5xl0'3V) 

REAL 

I  MAG 

REAL 

I  MAG 

S06 

03 

cc 

95 

-60.4619 

-1?. 42679 

M2. 94669 

♦29.99758 

.19806356 

-  .424870 

♦1.508699 

♦  .200000 

T08 

9Tb 

107 

-??. 01400 

-17.24408 

M5. 22787 

♦20. 34345 

.2006144 

.6929016 

♦  .4037744 

♦  .6932108 

uio 

* '  B 

115 

-  - 

-77.43989 

-20.49544 

M9. 73739 

^29. 20841 

.9140537 

-  .6743707 

♦  .5295426 

*  .4071000 

VI? 

flTB 

101 

-  6.23398 

-115.4021 

. 

♦38.51 152 

i?9. 45295 

.2606750 

-  .294674 

♦  .3271609 

♦2.649029 

y  14 

4TB 

no 

-24.31728 

-  6.543068 

♦56.38387 

♦29.34473 

.06778622 

.4046351 

♦  .2581218 

♦  .05179789 

X 16 

_ 

?tb 

11? 

-58.65565 

-44.09909 

♦67.08458 

♦37.64018 

.3214495 

-  .6434123 

♦  .2581229 

♦  .4217677 

1 1 BS  *  Interval  between  samples  used  by  Prony's  method 
for  calculating  the  poles  (Tg-200/255  nsec). 

2 

ITS  *  Starting  point  of  the  window  used  by  Prony's 
method  for  calculating  the  poles. 

The  first  pole  p?ir  rcrrnspon<*s  to  the  antenna  and 
the  second  to  the  tunnel. 


84 


C.  Results 

“  \ 

The  tunnel  resonance  as  given  by  Prony's  method  (20-30  MHz) 
corresponds  very  much  to  the  calculated  one  of  25  MHz.  The  deviation 
of  the  measured  resonance  at  different  recording  positions  is  due 
to  the  changing  height  of  the  tunnel  along  its  path. 

looking  at  Figures  19  and  20  we  observe  that  the  first  zero 
crossinq  of  the  reconstructed  tunnel  response  after  its  deviation 
from  the  filtered  waveform  occurs  at  about  90  nsec.  This  corres¬ 
ponds  to  the  arrival  time  of  the  tunnel  response  and  compares  satis¬ 
factorily  to  the  calculated  one  of  100  nsec. 

In  order  to  determine  the  relative  position  of  the  tunnel, 
grey- level  mapping  (see  Reference  '??])  of  the  reconstructed  waveforms 
for  each  traverse  could  be  performed.  This  mapping  technique  places 
the  waveforms  of  the  traverse  in  a  two  dimensional  coordinate  system 
of  depth  (time)  versus  waveform  recording  position.  Each  coordinate 
point  corresponds  to  a  waveform  voltage.  Voltages  between  successive 
positions  are  interpolated  by  using  Lagrange  polynomial  approximations 
(22).  Thus,  a  continuous  two  dimentional  array  of  voltages  is  obtained. 
A  computer  routine  written  by  Stapp  translates  the  array  voltages 
into  characters  of  quantized  darkness  levels.  The  highest  voltage 
corresponds  to  the  darkest  character  (black),  and  the  lowest  to 
white.  These  characters  could  then  be  plotted  on  a  line  printer 
to  obtain  a  grey- level  map  of  depth  versus  waveform  recording  position. 

The  grey- level  map  indicates  the  arrival  time  of  the  tunnel 
response  at  each  echo  recording  position.  As  we  approach  the  tunnel, 
the  arrival  time  of  its  response  becomes  shorter.  Such  behavior 
creates  a  hyperbola  on  the  map  with  its  peak  corresponding  to  the 
position  with  the  shortest  distance  from  the  tunnel  (for  a  leveled 
surface  this  position  is  right  over  the  tunnel). 

Figures  38-40  show  generated  maps  for  the  lower  traverse  of 
the  plan  view  given  in  Figure  37.  As  seen,  multiple  hyperbolas 
are  created.  Rut.  these  are  ficticious  due  to  the  reconstruction 
process.  The  hyperbola  corresponding  to  the  tunnel  is  the  one  be¬ 
ginning  at  about  90  nsec  (in  the  middle  of  the  map).  This  is  deter¬ 
mined  from  the  first  zero  crossing  of  the  reconstructed  waveform 
after  its  deviation  fran  the  filtered.  Figure  38  indicates  the 
map  with  interpolation  in  the  position  axis  performed  after  recti¬ 
fication  (Fold  First).  This  rectification  process  was  found  necessary 
in  Reference  (2?|  to  clearly  depict  the  tunnels  in  the  grey- level 
mapping.  It  is  not  needed  after  using  the  processing  techniques 
outlined  here.  The  black  level  was  set  at  3.5x5mV.  Other  descriptive 
material  on  these  figures,  not  needed  in  the  present  discussion, 
is  retained  for  the  deeply  interested  reader  and  follows  the 
definitions  given  in  (22).  In  Figures  39  and  40  interpolation  was 
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Figure  38.  Happing  of  the  lower  traverse  over  the  tunnel 
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F  igure  39.  Mapping  of  the  lower  traverse  over  the  tunnel 
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i iqure  40.  Mappinq  of  the  lower  traverse  over  the  tunnel. 
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performed  before  rectification  (Fold  Last)  and  with  the  black  level 
set  at  3.5x5  and  5x5mV,  respectively.  The  actual  tunnel  position 
at  recording  coordinate  RIO  (position  over  the  tynnel)  is  shown 
in  all  maps  for  comparison.  This  is  exactly  the  position  shown 
by  the  maps  at  this  coordinate.  But,  the  hyperbola  does  not  peak 
at  RIO  as  it  would  be  expected.  It  peaks  at  P6-Q8,  although,  it 
does  not  differ  much  frgm  position  RIO.  This  distortion  can  be 
well  attributed  to  a  10°  ground  slope  not  accounted  for  in  the 
maps,  which  could  have  made  the  shortest  distance  from  the  tunnel 
to  be  at  position  P6. 

In  Figure  41  a  grey-level  map  for  the  top  traverse  is  shown. 

The  actual  tunnel  position  at  coordinate  Ulo  is  the  same  as  the 
one  given  by  the  map.  Again,  as  discussed  previously  the  hyperbola 
peaks  at  positions  S6-T8. 

It  should  be  noted  that  in  all  maps  the  height  of  the  tunnel 
at  each  recording  position  corresponds  to  the  width  of  the  hyperbola. 

Our  results  as  compared  to  measured  ones  are  very  good. 

In  some  instances  measured  and  calculated  results  have  a  perfect 
correspondence,  while  in  others,  they  are  very  close.  The  important 
point  to  be  understood  here  is  that  our  process  is  indeed  capable 
of  identifying  the  tunnel.  A  comparison  o*  our  maps  with  those 
generated  by  Stapp  enhances  the  successfulness  of  the  process. 

Figures  4?  and  43  present  an  example  of  the  maps  given  by  Stapp 
for  the  top  and  lower  traverse  in  Figure  37,  respectively.  They 
were  produced  by  direct  grey- level  mapping  of  the  echoes.  As  seen, 
the  expected  hyperbola  w'ich  identifies  the  tunnel  is  not  clearly 
indicated.  More  likely,  the  maps  give  a  confused  picture  of  the 
tunnel  and  limited  conclusions  can  be  made. 

Stapp  was  forced  to  use  various  data  processing  techniques 
such  as  taking  the  absolute  value,  average  nearest  neighbor  waveforms 
and  adjusting  the  level  at  which  the  waveforms  were  clipped  etc. 

Even  then  the  results  while  acceptable  left  much  room  for  improvement. 
However,  by  including  the  physical  mechanisms,  i.e.  natural  resonances, 
in  the  process,  these  rather  arbitrary  steps  were  not  necessary 
and  vastly  improved  maps  were  obtained. 
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Figure  4?.  Mapping  of  the  lower  traverse  over  the  tunnel 
as  given  by  Stapp  (22). 
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CHAPTER  VI 

SUGARY  AND  CONCLUSIONS 


The  tunnel  identification  process  consisted  of  three  main 

tasks: 


1.  Recording  of  the  echoes  and  calculation  of  their  poles 
by  Prony' s  method. 

2.  Processing  of  the  echoes: 

a)  Removal  of  the  antenna  pole  pair  and  other  undesired 
resonances  (Chapter  II). 

b)  Reconstruct  ion  of  the  tunnel  response  (Chapter  III). 

3.  Determination  of  tunnel's  structure,  depth,  and  relative 
position  (mapping). 

He  concentrated  on  the  processing  of  the  echoes.  The  processing 
procedure  was  thoroughly  analyzed  and  tested  on  theoretical  (Chapter 
IV)  and  measured  (Chapter  V)  tunnel  echoes.  Methods  for  identifying 
and  determining  the  tunnel's  structure  and  depth  were  Outlined. 

The  reconstruction  process  was  used  for  combating  the  clutter  problem. 
The  relative  tunnel  position  could  be  obtained  by  employing  the 
grey- level  mapping  technique. 

Our  results  from  the  measured  echoes  were  very  good.  The 
tunnel  position,  depth  and  height  were  dearly  indicated  in  the 
constructed  maps,  and  they  did  correspond  to  the  actual  ones. 

Thus,  it  was  demonstrated  that  our  processing  can  indeed 
identify  and  provide  us  with  concise  information  about  an  under¬ 
ground  tunnel  from  its  backscattered  responses. 

The  processing  capability  was  very  accurate  as  long  as  the 
calculated  waveform  poles  by  Prony’s  method  are  within  reasonable 
accuracy  (especially  the  imaginary  part).  The  inaccuracy  of  Prony's 
method  when  dealing  with  data  dominated  by  clutter  and  noise  limits 
our  processing  capability.  Once  the  proper  waveform  poles  are  known, 
the  identification  process  can  be  performed  within  a  great  accuracy. 
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An  important  part  of  the  processing,  which  can  be  performed 
in  real  time,  is  that  the  structure  and  depth  of  the  tunnel  can 
be  determined  before  mapping.  Mapping  can  give  us  a  better  overview 
of  the  relative  tunnel  position.  But,  It  requires  extensive  computer 
t ime  and  storage,  which  is  undesired.  Furthermore,  the  construction 
of  maps  needs  many  recorded  waveforms.  This  can  be  a  problem  since 
the  ground  surface  may  not  be  smooth  enough  for  the  antenna  to  be 
properly  positioned  in  many  spots. 

It  should  be  noted  that  our  processing  of  the  echoes,  although 
it  was  used  for  identifying  tunnels,  is  very  general.  It  can  be 
applied  to  any  target  backscattered  waveform  for  isolating  the  target 
response  and  combating  clutter  and  noise. 
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APPENDIX  A 

DERIVATION  OF  A  DIFFERENTIAL  EQUATION  FOR 
EXTRACTING  ONE  COMPLEX  CONJUGATE  POLE  PAIR 


Let  us  assume  an  original  waveform,  r(t),  with  a  Laplace  repre 
sent  at  ion  of 


R(s)  »^[r(t)]  . 

In  general  we  can  express  R(s)  as  follows: 

R(S)  ,  - -  (Al) 

(s-s,)(s-Sj)D(s) 

♦ 

where  (S|,Sj)  is  the  pole  pair  to  be  removed. 

Multiplying  both^sides  of  Equation  (Al)  by  the  representat ion 
of  the  pole  pair  (s,,Sj).  we  have 

Rp(s)  -  (s?-?sRe(s1)*|s1|?)R(s)  «  .  (A 2) 

R_(s)  is  the  Laplace  representation  of  the  filtered  waveform.  Its 
poles  are  described  by  D{s). 

Transforming  Rp(s)  to  time  domain  we  obtain 
rp'*)  -(j£  -  »«<*,)  ^  *  |s,l)r(U  (»3) 

where  rp(t)  *^_1[Rp(s)]  . 

Equation  (A3)  indicates  the  analog  implementation  required 
for  extracting  the  complex  conjugate  pole  pair  (s.,s.).  It  is  equiva 
lent  to  the  digital  one  given  in  Equation  (14)  (see  Chapter 
II). 

In  a  similar  manner  we  can  derive  differential  equations  equiva 
lent  to  the  difference  equations  given  in  Chapter  II  and  III. 
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APPENDIX  B 

MAIN  FORTRAN  PROGRAM  FOR  CALLING  THE  POLE  EXTRACTION  AND 
RECONSTRUCT  ION  ROUTINES,  MANIPULATING  DATA  AND  PLOTTING 
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APPENDIX  C 

THE  POLE  EXTRACTION  SUBROUTINE 
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APPENDIX  D 

THE  POIE  RECONSTRUCTION  SUBROUTINE 
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APPENDIX  E 

THE  TRANSMISSION  LINE  COMPUTER  MODEL  WITH 
ADDITIONAL  RESULTS 
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( •  rHL'’l«rHL;T*<I.K*ft) 

t<  FrtfcUAaMtEOA*  I  I  .  F'RY  ) 

lii  t|li|*blkT*H.rl'l>l 

tn  SloA-SllAK I .n'ft> 

•  V  '.'u* -x«l -F*l  .•  T|  •  I  .  t-W  ) 

i.  Ff  s»t-w*c(  i/« ;g«pi  )  1*1  .E-vi 

.  I  •  Ft  |  MV)  •Ml,  EPS 

...  .IV  FOUilAfl  ll>2U2«).A> 

,OA|oAA*^*PI«F«kCA 
i»  (jWfcii  AT*.1  *P  I  *rn  hC7 

iV  v  »»  r..»;-J',AL  II  ‘E  i>|NOCM  I '  ES I  RET)  (VSEC) 

•  <  UX  »»  .  •  ■  i  .  <  r  r JED.  POINTS  I*  THE  TRANSFER  FUNCTION 

LoC  »>l-r  •f-Rh'HJfc'  CY  |TCHEl'E»T  LEVF1. 

It  NPMS-JSe 

Tv  r:p  s*M2 

o.  nzp*ij 

ol  Nr  E  to*  4k  v® 

t.  v  •  •  •  •  • 

t»  V.  •  •  oTfcR  1  IN  ATI  ON  OF  TI'*E  ATO  FftFOUFNCY  INTERVAL 

•  •  •  •  • 

ct  ■  ..!•>' EG**  •Trklo 

■  N4»  .  irliEO/4 

!«  T'.i"i»T ■)/  II'. 

■  v  i.r  •(  )/..  )*I.Ff*3 

v»  I0*|4i/lf  PS-1  .» 

vl  r  »|  •(firk£0/T!  >•  1  .EPV 

v.  >l|  •  r«  I .  IMc 

Vj  nA.;*i>»ih  A 

V*  C  »HA*I  "TETNA  Sl'IKT  kFSI  STOP 

v.  sjiiwt  resistor 

Vo  I  r  (  N 1 0 A .  ) .  )  GO  TO  *<« 

f.A*Ai.t(»«A/(2*sirj  m 
V  .1  4  mi  •  l  . 

w  I F < S I l/T . EG . i* .  1  GO  TC  flj 

U  ..•A.‘.MRST/«2*Stcn  J> 

III  lj  Hl»  •  I  • 

li.  RS*SlO  IA*(  MV/FP?  ) 

Ik  J  V.  •  •  •  •  • 

Ik  »  c  •  •  :  ALCI  ’I.  AT  I  OH  OF  ARRAY  VOUT 

Ut>  v.  •  •  •  •  • 

i.o  JO  I.V  I-Z.NFkEO 

i.  ■  .  •DEIFRt'I  *  E  ANTED*  a  TARGET  I  fPEHAr  CF  (Z a  A  ZT  ) 

I.  •  CARL  *k.A  Pl.A  ( k’.  .OHEGa  ) 

Ik  V  POLA«;vPLA(S|r,A.  flEr.AA  ) 

III.  CFOLA-CCNJGC  nn.A  > 

ZA«(  nSA*CFkFO)/<  a>uFO-POtA  )*(CF«FO-CTO|.A>  > 

II.  POLT^J'TLX  (S IGT  ,  O  COAT  I 

Mj  CPOL  KUIJGI  POLT  I 

114  Z  i  •(  mS7*CFRF(.) /(  (CFRFO-POI  7  )«(OFRFO-Cr,CLT  I ) 

llj  o  •  •  T  A  m  F  i  IfcPF»Am.f  AT  INPUT  TFPi'I  Tit  S  (ZTFr*' ) 

II.,  o  >»2t»I.IFE  CHARACTERISTIC  I XPFGAHCF 
117  o  »»(,A>A»I  ROPArAlll*’  COST  ATT 

Ilk  ZARG*UTfcrA*RU 

llv  2 AUG' i*  fl:GA*EPS 

l.i  Zk'Anl  ..  I  P|.A(RG,ZAIir.)A:*'P|.X(Sir.l,A./ARrj  » 
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14  1  4.)«A'SM<'»  UUAl*. ) 

I44  OA«A<  :PLX(»f  ./AWUl/2.1 

•  4-  hOAJA-Kfc  At.  CC/HA) 

I.'-  lAjAV.  A.  A  i  .IAO(  <  A  A» 

i.  a  •  •  tliil  Hh  alCULaTIOn  iNl.>  •ir-t-LFCTIDH  a 

l^(.  A  UANSIJSMaM  coErr  let  hits 

147  IHIECHt.EO.t  )  on  TO  K 

i4>t  *ts*)77. 

I<V  ZA.’»<ZA»«:,M/(/A»/.)) 

I..0  «;*L«iAkl/</A.,»>ltS) 

IJI  ZnfcSA*UA*kEJ J/CZA.rfS) 

1.4  ziK»«(^-ir^ir>A)/i  /i.t-SA*.:.’) 
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I  JO  C  l*K*l  • 
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15  I  105  MM*  I  . 

154  a  .  .  iiAi'fi.  fo'icti  • 

•  5  j  «u.i  ( I  >«<mi  i  j/ c i  j) 
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2  >  f»«M  f  I*  "»,(  IX.'Alt*  I  t  ’.1(1)1  ') 

2  7o  .’Af 

..  If  U  tAl.L  Of*  •  C  I  AF.EOI  )  V  *.«:  1  "4 
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2a i  »..-4  coniitnie 
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2b  V  CALL  411  All  I  rL  I,M) 

2*4  o  •  •  •  •  * 
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2V*  0  •  •  •  •  • 

2*4  PAUSE 
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2*o  .  14  It.f  (l  (  I  ) 
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.  *  A  ,15  '..If  o  (  I  .f  (.I  |  I 
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-H  I  "Cf  I  •*%.»'  |«l,2 

•»  .  all  ’Li  4  <  4.b. I .  *?b.-J> 

Jl>.  ALL  A  1*1...  ..Ill 

J*’*.  C2.,I.,1,J 

J»>b  IMl.fcCJ..?)  C.C  ro  212 

C ALL  «XIS<e.  I MTIMC  (H8EC>.*ll,ll>.*M).«C..1Kf, 

jt-  *  *.>..*'> 

J*t  CALL  LI!. I  ir.T  .«>..^..I!JF<;..,..Tf:,.,,P5..*.(*> 

■>»  ■»  CO  i  J  J  I  j 

J  1  *  .12  1 .  h  •  I  . 

Ji»  CALL  A«irc.>.,  . .  1 1  IfTI  ME  (  •rrC) 1 1  .11*.  .V  *.  .^TP'.T^n, 

J  I  2  \l..4*> 
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J  I  A  .  I  J  CO?’.  ’I.ofc 
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J»e  IU  •  •  !  ••  1. 11.:  Si;  r. I*.) 

J  I  4  .1C  »it<  •  I  . 

J'b  v  ALL  iYMi/L«2..A...II.lfcL*“.v»A..  ,*  > 
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->.  •  CALL  MaJ  St"  <2.  ,©.rt,.  1 1  ,v,*.  ,J 
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J*./  ALL  :  i!  fcxi  2.a. -j.2*)  ..  1 1  .1 1  .w. .  1 1 
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ALL  -.i!  M2.6,  •>.  7,  .  1  1 ,  siC'A.oi*.,  |p?> 

JJi  ALL  SYl  l-i  L«2.».4.v. .  1 1 ,  irtSTA.oit.  ,|A) 
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ALL  <oi Lc  fct  ...  .  4  . .  I  1 ,  :!.•*...?) 
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Jbb  Oft*  ilt!i  SSl*.4> 

Abo  COA.PILX  IJATA  l*:> 
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Additional  Results  on  the  Transmission  Line  Model 


We  present  here  the  results  of  five  more  sets  of  parameters 
used  in  the  transmission  line  model  for  calculating  the  target 
(tunnel)  pole  pair  and  its  depth.  The  first  echo  responses  are 
given  in  Figures  El,  E3,  E5,  E7,  E9  and  their  respective  filtered 
waveforms  (double  antenna  pole  pair  extracted)  and  reconstructed 
target  responses  in  Figures  E2,  E4,  E6,  E8  and  E TO .  Attention  should 
be  given  to  the  first  set  of  parameters.  In  this  case  Equation 
(61)  is  not  satisfied  for  the  target  response.  Thus,  the  imaginary 
part  of  the  calculated  target  pole  (Table  El)  deviates  from  the 
actual  one  (Table  2).  This  situation  becomes  more  cirtical  as  the 
target  resistance  Increases  with  respect  to  Z  .  The  same  is  also 
true  for  the  antenna  pole  pair.  Under  such  conditions  we  can  still 
calculate  the  target's  depth,  but  not  its  structure.  It  should 
be  noted,  though,  that  given  the  transmission  line  parameters  (Z 
and  R  )  we  could  finally  calculate  the  proper  target  resonances  To 
determine  the  target's  structure.  All  sets  have  the  following  common 
parameters: 

Conductivity  of  o  *  5.x10*^U/m. 

Line  characteristic  impedance  of  ZQ  ■  377. ♦jO. 

Generator  internal  resistance  of  R$  a  377. Q 

Relative  permitlvity  of  cr  *  1. 

Relative  permeability  of  pr  ■  1. 

Line  length  (target  depth)  of  d  a  30  m. 


Table  El 

MODEL  PARAMETERS 


WVFRM 

TARGET* 

P0LE(x106) 

ANTENNA* 

POL  E ( xlO6) 

ra 

(ohms) 

rt 

(ohms) 

tb 

(nsec) 

GAUSSIAN 

PULSE  WIDTH 

Fig.  IE 

-50*j40 

-50* J50 

100 

600 

.39064 

6  nsec 

Fig.  ?E 

-30* J60 

-40* J70 

312 

416.67 

.61043 

6  nsec 

Fiq.  4E 

-60*j80 

-50* J90 

250 

208.33 

.61043 

6  nsec 

Fiq.  6E 

-40* J80 

-90* j 100 

277.78 

625.0 

.73251 

6  nsec 

Fig.  8E 

-8CHJ100 

- 100+J 120 

250 

312.5 

.48834 

3  nsec 

*  The  real  and  imaginary  parts  of  the  poles  are  in 
Nepers/s»c  and  H *,  respectively. 
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The  real  and  Imaginary  parts  of  the  poles  are  given 
in  Nepers/sec  and  Hz,  respectively.  The  first  two 
pole  pair:  correspond  to  the  double  antenna  pole 
pair,  and  the  third  to  the  target. 
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Figure  El.  First  echo  response. 
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Figure  E 2.  Extraction  of  antenna  double  pole  pair  and  reconstruction 
of  target  response  of  the  signal  in  Figure  El. 


gure  E4.  Extraction  of  antenna  double  pole  pair  and  reconstruction 
of  target  response  of  the  signal  in  Figure  E3. 
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Figure  E6.  Extraction  of  antenna  double  pole  pair  and  reconstruction 
of  target  response  of  the  signal  in  Figure  E5. 
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Figure  E8.  Extraction  of  antenna  double  pole  pair  and  reconstruction 
of  target  response  of  the  signal  in  Figure  E7. 


Figure  E10.  Extraction  of  antenna  double  pole  pair  and  reconstruct 
of  target  response  of  the  signal  in  Figure  E9. 


